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1.  INTRODUCTION 


The  Plasma  Physics  Division  of  the  Naval  Research  Laboratory  has  long 
played  a  leading  role  in  advancing  the  physical  understanding  of  intense 
charged  particle  beam  generation  and  acceleration.  A  particularly 
noteworthly  aspect  of  NRL's  contribution  to  the  field  is  presented  by  the 
series  of  Innovative  techniques  proposed  by  the  Lab  to  accomplish  the 
acceleration  of  multi -kilo amp  electron  beams  to  energies  well  above  100 
MeV.  Among  these  suggested  schemes  is  the  Collective  Particle  Accelerator 
(CPA)  concept  Introduced  by  Or.  Moshe  Friedman  (Reference  1)  in  1977  and 
experimentally  Investigated  by  him  for  about  three  years  thereafter. 
Although  initial  results  were  negative,  recent  theoretical  analysis  by 
Dr.  Peter  Palmadesso  (Reference  2)  indicates  that  some  simple  changes  in 
the  device's  operating  parameters  could  produce  a  positive  result.  Given 
the  potential  significance  of  such  an  accelerator  to  various  programs 
underway  In  the  Navy,  it  is  highly  desirable  to  resolve  the  question  of 
the  CPA's  feasibility.  Reconstructing  the  now  dismantled  experiment  to 
test  the  suggested  theoretical  modifications  would  be  an  extremely  costly 
and  time-consuming  undertaking  which  would  be  better  justified  if  an 
Independent  verification  of  the  theorized  Improvements  were  available. 
For  this  purpose.  Mission  Research  Corporation  (MRC)  conducted  two- 
dimensional,  relativistic,  electromagnetic,  particle-in-cell  simulations 
using  a  model  of  the  CPA  with  realistic  operating  parameters.  The  results 
of  these  simulations  have  pinpointed  the  probable  cause  for  the  previous 
problems  and  indicate  possible  solutions  which  could  spell  success  for  a 
renewed  experimental  effort. 

In  the  technical  discussion  that  follows,  the  general  research  topic 
Is  described  and  the  specific  objective  of  this  work  is  delineated.  The 
discussion  begins  with  a  brief  review  of  the  experimental  CPA  device  in 
question.  The  underlying,  general  theory  is  spelled  out  and  the 
fundamental  operative  equations  are  listed.  Thereafter  follows  a  detailed 


picture  of  the  actual  experimental  apparatus  constructed  at  URL  to  test 
the  feasibility  of  the  concept.  The  observed  negative  results  are  then 
described.  At  that  point,  the  specifics  of  this  effort  are  initiated  with 
a  reduction  of  the  actual  physical  device  to  a  numerical  model  which  is 
both  physically  valid  and  computationally  tractable.  Particular  attention 
is  paid  to  the  assumptions  which  are  made  to  arrive  at  the  model.  These 
assumptions  determine  which  phenomena  can  or  cannot  manifest  themselves  in 
the  simulations.  Having  fixed  the  subject  of  the  study,  the  discussion 
then  turns  in  Section  4  to  the  computational  tool  which  was  applied  to  its 
analysis.  For  this  work  MRC  used  its  2-1/2-D,  electromagnetic, 
relativistic  MAGIC  code.  A  complete  description  is  given  of  the  code's 
capabilities  and  limitations.  Section  5  presents  the  detailed  results  of 
the  five  major  sets  of  simulation  runs  conducted  with  MAGIC  to  diagnose 
the  CPA  device.  Numerous  figures  are  included  to  provide  as  complete  as 
possible  a  picture  of  the  various  physical  phenomena  at  work.  The  major 
experimental  observations  are  duplicated  and  their  causes  are  carefully 
diagnosed.  Clear  conclusions  are  drawn  in  Section  6. 

A  short  appendix  is  included  at  the  end  of  this  report  to  document  a 
concise  but  important  piece  of  software  developed  by  MRC  for  the  NRL  VAX 
during  the  course  of  this  work.  This  software  adapts  the  NCAR  Graphics 
System  hosted  on  that  machine  to  respond  to  commands  for  the  much  more 
widely  used  DISSPLA  Graphics  System.  This  feature  has  now  been  made 
available  to  all  users  of  the  machine. 


2.  THE  EXPERIMENTAL  DEVICE 


Mission  Research  Corporation  has  applied  its  MAGIC  simulation  code  to 
examine  specific,  critical  physics  questions  involving  the  technical 
feasibility  of  NRL's  Collective  Particle  Accelerator  (CPA).  These 
simulations  now  provide  critical  additional  information  to  the  parties 
concerned  at  NRL  to  determine  the  possible  future  of  a  revived  CPA 
experimental  effort.  If  the  CPA  is  again  brought  to  life,  MRC  proposes  to 
continue  its  computational  support  of  that  program  to  optimize  its 
performance.  This  section  of  the  report  reviews  the  NRL-CPA  experiment 
and  its  chief  results. 

An  approximate  drawing  of  the  CPA  device  experimentally  tested  at  NRL 
by  Dr.  Moshe  Friedman  appears  in  Figure  1.  The  essential  elements 
consisted  of  two  opposing  e-beam,  foil less  diodes  installed  at  opposite 
ends  of  a  five-meter-long,  4.76  cm  diameter  metal  drift  tube.  Both  diodes 
were  fed  from  the  same  1.0  megavolt,  100  nanosecond  pulsed  power 
generator.  However,  a  transformer  Installed  in  the  feed-line  branch 
leading  to  the  hollow  cathode  diode  introduced  an  impedance  mismatch  in 
that  line  which  resulted  in  an  effective  diode  voltage  at  that  end  of  1.5 
MV  compared  with  the  1.0  MV  of  the  opposing  rod  diode.  The  hollow-beam 
diode  consisted  of  a  4.0  cm  diameter  tapered  graphite  cathode  inserted 
into  the  flared  anode  structure  at  one  end  of  the  drift  tuber  The  anode 
diameter  Increased  from  4.76  to  15.24  cm  over  an  axial  distance  of  12.7 
cm.  This  flaring  was  accomplished  in  two  parts:  first,  a  diameter 
Increase  from  4.76  to  7.0  cm  over  a  5.1  cm  axial  distance,  and  then,  an 
Increase  from  7.0  to  15.24  cm  over  the  remaining  7.6  cm.  This  two-staged 
anode  flaring  was  exactly  repeated  on  the  opposite  end  of  the  drift  tube 
where  the  rod  diode  was  placed.  The  axial  positioning  of  the  cathodes  was 
not  a  fixed  parameter  in  the  experiments.  The  electron  beam  produced  by 
this  cathode  had  an  outside  diameter  of  about  4.0  cm  and  a  radial  beam 
thickness  of  about  0.2  cm.  It  carried  a  peak  average  current  of  about  20 


kiloamperes.  The  tip  of  the  rod  cathode  typically  projected  one  or  two 
centimeters  into  the  drift  tube.  A  graphite  ring  less  than  1.0  cm  in 
cross-section  hugged  the  inner  drift  tube  wall  in  the  vicinity  of  the  rod 
tip  to  act  as  a  beam  stop  for  the  hollow  electron  beam.  The  axial 
position  of  this  graphite  beam- stop  was  again  not  a  fixed  experimental 
parameter.  In  some  cases  it  was  placed  a  short  distance  in  front  of  in 
some  cases  a  short  distance  behind  the  rod  cathode  tip.  The  rod  cathode 
emitted  a  pencil  e-beam  of  diameter  0.5  cm  with  a  peak  current  of  2  kA  and 
a  current  and  voltage  rise  time  of  about  5.0  ns.  In  the  actual  CPA,  the 
hollow  electron  beam  would  experience  a  density  modulation  sinusoidally  to 
zero  at  a  predetermined  frequency.  Uniformly  positioned  several 
centimeters  apart  down  the  entire  external  length  of  the  drift  tube  were  a 
series  of  magnetic  field  coils  which  could  create  either  a  constant  15 
kilogauss  Bz  or  an  axially-rippled  B- field  consisting  of  an  average  22.5 

2itZ 

kilogauss  component  to  which  is  added  a  ripple  of  the  form  Bx  .cos  (  „  ) 

where  Bj  -  7.5  kG  and  L  =  60.0  cm.  (Note  that  the  uniform  field  coil 
spacing  was  compressed  for  the  last  several  coils  at  both  ends  of  the  tube 
to  provide  intensified  Bz-fields  behind  the  respective  diodes.) 

To  describe  the  physics  of  the  CPA  device,  we  rely  on  excerpts  from 
Reference  1.  In  that  report,  it  is  pointed  out  that  the  axial  electric 
field  component  at  r=0,  produced  by  an  annular,  unneutralized, 
magnetically-focused  intense  relativistic  electron  beam  { IREB)  propagating 
through  a  drift  tube  of  radius  R  can  be  approximated  by 
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where  Q  is  the  beam's  charge  per  unit  length  (Q  =  I/v^).  I  is  the  beam 
current,  rj,  is  the  beam  radius,  and  v  is  the  electron  velocity.  If  the 
drift  tube  walls  are  covered  with  a  thin  dielectric  layer  of  thickness, 
6R,  and  permeability,  e,  one  finds 


2  ir  3t  rb  2ttB2  e-1  R  at  2ire0  rfa  az 


Furthermore,  by  choosing 
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Equation  (2)  reduces  to 
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Now  assume  a  charge  density  modulation  of  the  IREB  of  frequency  and 
wavelength,  f  and  \  respectively,  such  that  xf  =  v  and 
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If  an  axially,  sinusoidally  rippled  magnetic  field  with  ripple  wavelength, 
L,  is  now  imposed  in  the  drift  tube,  the  annular  beam  will  follow  the 
B-field  lines  and  have  an  axially  varying  radius  of 


-  r  +  r,  cos  — —  :  r,  <  r_ 


(6) 


where  r0  is  the  mean  IRE8  radius  and  rl  is  the  amplitude  of  the  radius 
ripples.  Inserting  {5)  and  (6)  into  (4)  along  with  the  assumption  that 
the  beam  velocity,  v,  is  approximately  equal  to  the  speed  of  light,  c, 
yields 


Thus,  there  are  two  "waves"  with  phase  velocities 


v  =  v -  ,  for  the  forward  wave,  and 

■  X  +  L 

('8) 


V  o  =  -v -  ,  for  the  backward  wave. 

X  -  L 


The  amplitude  of  these  "waves"  is  given  by 


Since  only  the  backward  wave  can  have  a  phase  velocity  approaching  c,  only 
the  backward  wave  can  accelerate  electrons.  Note  that  both  waves  can 
accelerate  ions  but  ion  acceleration  is  beyond  the  scope  of  this  proposed 
research. 


In  the  experiments  conducted  to  date  with  the  device  in  Figure  1,  the 
hollow  beam  density  was  not  modulated  so  acceleration  could  not  be 
tested.  However,  for  the  beam  and  rippled  Bz  that  was  in  place,  L  * 
0.60  meter,  rx/r0  »  0.17,  and  Q  *  I/c  =  2xl0V3xl08  *  6.67xl0-5 
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coulomb/meter.  Assuming  that  the  hollow  beam  could  be  density  modulated 
with  x  -  2L,  then  Equation  (8)  predicts  a  phase  velocity  of  c  for  the 
backward  wave  which,  with  this  value  for  Q  inserted  into  Equation  (9) 
yields  530  kV/meter  as  the  effective  axial  acceleration  field.  Over  the 
full  5  meter  length,  this  amounts  to  a  2.5  MV  energy  gain  which  should  be 
easily  observable  experimentally.  It  is  worth  noting  as  well  that  an  IREB 
propagating  through  a  rippled  magnetic  field  can  produce  significant 
microwave  radiation  if  the  field  strength  is  below  a  critical  value  given 
by 

B  .  il  —  Y  •  (10) 

c  L  e 

For  the  experimental  parameters  of  L  *  0.6  and  Y  =  1  +  1.5/. 511  =  3.9, 
Bc  becomes  about  0.7  kG  which  is  far  below  the  average  15  kG  baseline 
axial  guide  field.  It  is  also  important  to  note  the  axial  and  temporal 
scale  lengths  in  the  experiment.  The  relativistic  electron  gyroradius  is 
given  by 

rQ  =  1.7xl03  (Y2  -  1)1/2  B”1  cm  (11) 

so  that  in  the  peak  axial  field  of  30  kG,  rg  -  0.22  cm.  Similarly,  the 
relativistic  electron  gyroperiod  is  found  to  be  about  5x10" 11  second. 

In  the  actual  experiments  that  were  run  at  NRL,  the  hollow  beam  was 
not  density  modulated  and  the  central  pencil  beam  apparently  did  not 
propagate  down  the  full  length  of  the  tube  (Reference  3).  Thus,  the 
experiment  never  proceeded  to  the  point  of  testing  the  acceleration 
concept  itself.  Specifically,  with  the  20  kA  hollow  beam  in  steady  state 
propagation  down  the  full  length  of  the  five-meter  tube  (confirmed  by 
damage  patterns)  the  pencil  diode  voltage  was  ramped  up  from  0  to  1.0  MV 
over  about  five  nanoseconds.  The  pencil  diode  current  rose  likewise  from 
0  to  2.0  kA  over  about  the  same  five  nanoseconds.  The  pencil  beam  was 


never  observed  to  reach  the  other  end  of  the  tube  in  the  presence  of  the 
hollow  beam.  In  fact,  magnetic  field  probe  measurements  made  just  50  cm 
downstream  from  the  pencil  diode  failed  to  detect  the  presence  of  the 
pencil  beam.  Therefore,  the  most  significant  immediate  question  addressed 
by  these  numerical  simulations  of  the  device  was  the  simple  one  of 
survivability  of  the  central  beam.  This  report  documents  the  numerical 
duplication  of  all  of  the  above  experimental  observations.  In  addition, 
it  was  found  that  the  central  beam  does  in  fact  propagate  down  the  axis  of 
the  tube  but  that  its  current  is  drastically  reduced  due  to  the  formation 
of  an  effective  virtual  cathode  right  next  to  the  emitting  tip  of  the 
central  rod  cathode. 

Before  describing  these  results  in  Section  5,  the  following  two 
sections  first  detail  the  numerical  model  of  the  CPA  which  was  used  in 
this  effort  and  then  explains  the  major  features  of  the  MAGIC  simulation 
code. 


3.  THE  NUMERICAL  MODEL 


From  a  practical  computational  standpoint,  it  was  clear  that  a 
reasonable  simulation  must  be  modeled  upon  a  device  significantly  shorter 
than  the  one  depicted  in  Figure  1  and  one  with  some  simplifications  in  the 
two  beam  injection  regions  as  well.  To  arrive  at  that  conclusion, 
consider  the  spatial  and  temporal  resolution  mandated  by  the  physics  of 
the  problem.  As  stated  in  the  previous  section,  the  gyroradius  for  the 
1.5  MeV  electrons  in  the  30  kG  magnetic  field  is  about  0.22  cm  with  a 
gyroperiod  of  5.0x10" 11  second.  Even  more  stringent,  for  the  1.0  MeV 
pencil  beam  these  quantities  become  0.16  cm  and  3.6  x  10" 11  second 
respectively.  The  maximum  numerical  gridspacings  in  either  dimension  and 
the  timestep  must  be  small  enough  to  resolve  those  electron  motions. 
(Luckily  the  beam  currents  are  low  enough  so  that  the  maximum  azimuthal 
magnetic  field  strengths  at  the  outer  radius  of  the  inner  and  outer  beams 
are  only  1.25  and  3.1  kG  respectively  which  are  small  compared  to  the 
axial  field.)  Thus,  the  axial  B-field  strength  peaks  midway  between  the 
two  diodes  so  the  largest  grid  spacing  in  that  region  must  be  less  than  or 
equal  to  0.16  cm.  Actually,  in  the  radial  dimension  there  is  even  more 
severe  a  restriction,  namely  that  which  demands  adequate  spatial 
resolution  along  the  electron-emitting  cathode  surfaces.  The  central  red 
cathode  has  an  emitting  tip  only  0.25  cm  in  radius  while  the  hollow 
cathode  emits  from  a  strip  only  0.2  cm  wide.  Clearly,  a  radial  cell 
spacing  of  0.16  cm  is  grossly  Inadequate  to  resolve  any  radial  beam 
profile.  To  rectify  this,  the  rod  cathode  Is  given  a  radius  of  0.32  cm 
and  the  hollow  cathode  emission  region  Is  widened  to  0.24  cm.  In 
addition,  the  radial  cell  spacing  near  those  electrodes  is  shortened  to 
0.08  cm.  Even  with  those  new  parameters,  the  rod  beam  emitter  is  only 
four  cells  wide  and  the  hollow  beam  only  three.  Widening  the  beam  radii 
further  seems  unwise  If  reasonable  resemblance  to  the  experiment  is  to  be 
maintained.  On  the  other  hand,  making  aR  smaller  than  0.08  cm  seems 
unwise  from  a  numerical  standpoint  due  to  the  limitation  of  computer  data 


LHW.WSfeil  ’/V  VA^X*7**rl  '7  **  -.’.  •.',  -.*  *•_■  *-•  -.“  -'■.■’-.r.jr.-.ir.  /-.  .-.  .-.  .-.  t.  ;\  i-  — ;  =-.  v.  .-. 


storage  for  the  numerical  mesh  points  and  to  the  resultant  time  step 
shortening  dictated  by  the  Courant  condition.  Summarizing  the  spatial 
grid  considerations  thus  far,  it  would  seem  desirable  to  utilize  aR  and  az 
spacings  of  0.16  cm  except  in  the  beam  emission  regions  where  a  resolution 
of  0.08  is  chosen  in  the  radial  dimension  to  provide  an  adequate  numerical 
beam  width  and  in  the  axial  dimension  to  ensure  accurate  determination  of 
the  normal  electric  field  for  use  in  the  Gauss'  Law  calculation  of  emitted 
charge.  With  regard  to  the  simulation  time  step,  it  is  desirable  to 
resolve  the  minimum  electron  gyroperiod  with  at  least  five  time  steps  to 
insure  physical  particle  dynamics.  This  permits  a  maximum  time  step  of 
about  7.2xl0"12  second.  The  Courant  numerical  stability  condition 
furthermore  demands  that  the  time  step.  At,  satisfies  the  expression 


1 

C2At2 


> 


J_  +  _I_ 

aR2  aZ2 


(12) 


For  the  smallest  cells,  near  the  cathode  tips,  aR  =  aZ  =  0.08  and,  thus, 
the  timestep  is  restricted  to  At  <  1.886  x  10~12  second.  Given  this 
temporal  constraint,  a  choice  of  At  *  1.40  x  10"12  second  was  deemed  quite 
reasonable. 


Now  consider  the  Implications  those  values  of  aR,  aZ.  and  At  have  on 
conducting  a  simulation  of  the  full  five-meter- long  CPA  device.  The  least 
difficult  starting  point  Is  the  radial  dimension.  For  convenience,  the 
grid  can  be  divided  into  five  regions  in  R  as  follows: 


(a)  From  R  *  0  to  R  *  0.48  cm,  use  a  uniform  radial  spacing  with  six 
cells  of  aR  *  0.08  cm  since  this  must  resolve  the  electron  motions  of  the 
central  beam. 


(b)  From  R  =  0.48  to  R  *  1.18  cm,  gradually  expand  the  spacing  in 
this  particle-free  region,  spanning  it  with  five  cells  whose  aR's  are 
0.10,  0.12,  0.14,  0.16,  and  0.18  cm  respectively. 
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(c)  From  R  =  1.18  to  R  =  1.70  cm,  as  one  approaches  the  domain  of  the 

outer  beam,  the  radial  cell  size  must  gradually  decrease  again.  This  is 

accomplished  with  four  cells  of  aR  =  0.16,  0.14,  0.12,  and  0.10  cm 

respectively. 

(d)  From  R  =  1.70  to  R  ==  2.10  cm,  again  use  a  uniform  radial  spacing 
of  aR  *  0.08  cm;  this  time  with  five  cells  to  span  the  beam  region. 

(e)  From  R  =  2.10  to  R  =  2.40  cm,  three  radial  cells  of  aR  *  0.09, 
0.10,  and  0.11  cm  fill  the  remainder  of  the  tube  Interior  out  to  the 
cylindrical  anode  wall. 

This  gives  a  total  of  24  cells  spanning  the  computational  region 
radially.  To  span  5.0  meters  in  the  axial  dimension  with  a  aZ  *  0.16  cm 
(ignoring  the  higher  resolution  needed  Immediately  near  the  cathode  tips) 
would  require  3125  cells  to  yield  a  total  of  22  x  3125  3  68,750 
two-dimensional  data  points  for  the  full  model.  A  simulation  of  that  size 
Is  doable  but  totally  Impractical.  Even  modeling  only  one  full  60  cm 
B-fleld  cycle  would  require  8250  cells  which  would  be  possible  but  not 
practical  if  more  that  just  one  or  two  simulations  are  desired.  To 

resolve  the  problem.  It  should  first  be  decided  to  examine  only  one 
magnetic  cycle,  L.  Furthermore,  L  should  be  shortened  from  the 
experimental  value  to  as  small  as  possible  a  size  without  sacrificing 
significant  physics.  In  essence,  this  means  that  L  must  remain  "many" 
times  the  value  of  the  tube  diameter  of  4.8  cm.  A  good  value  for  "many" 
Is  "5",  which  yields  an  L  of  about  24  cm.  Furthermore,  given  the 
experimental  uncertainties  of  the  ripple  field's  generation,  a  lowering  of 
the  peak  field  from  30.0  to  25.0  kG  hardly  seems  unreasonable.  This 
relaxes  the  minimum  cell  size  to  about  0.19  cm.  Such  a  value  for  aZ  would 
allow  the  use  of  130  cells  to  span  the  entire  CPA  model  axially.  The 
problem  has  now  assumed  practical  proportions  and  is  ready  for  detailed 
modeling. 


The  simulation  model  which  emerges  from  these  considerations  Is  drawn 
to  scale  In  Figure  2.  That  figure  depicts  blcw-ups  of  the  two  diode 
regions  at  either  end  of  the  tube.  The  treatment  of  the  radial  dimension 
Is  exactly  as  described  In  the  previous  paragraph.  The  axial  dimension  Is 

i 

a  bit  more  complicated.  Consider  first  the  hollow  cathode  at  the  left  end 
of  the  tube.  It  Is  no  longer  tapered  as  was  the  case  experimentally. 
Instead,  It  Is  a  2.02  cm  radius  cylinder  with  a  thickness  of  0.24  cm  and 
an  axial  depth  of  1.08  cm.  The  "backplane"  of  the  hollow  cathode  also 
defines  one  boundary  endplane  of  the  computational  region.  The  "open", 
freespace  region  on  this  endplane  between  R  =  2.02  and  R  *  2.40  cm 
represents  the  entrance  port  from  the  1.5  MV  power  transmission  line.  Its 
axial  recess  from  the  cathode  tip  of  more  than  twice  its  radial  width  is 
Important  for  ensuring  a  proper  numerical  representation  of  the  free 
boundary  potential  distribution.  The  volume  between  Z  =  0  and  Z  *  0.38  cm 
Is  spanned  by  two  axial  cells  of  width,  aZ  *  0.19  cm.  Over  the  final 
five  axial  steps  In  that  transmission  line  recess,  aZ  Is  decreased  from 
0.18  to  0.16  to  0.14  to  0.12  to  0.10  cm  until  the  end  of  the  cathode  shank 
is  reached  at  Z  *  1.08  cm.  An  axial  cell  width  of  no  more  than  0.08  cm 
(equal  to  the  best  radial  resolution)  is  desired  along  the  electron- 
emitting  cathode  tip  to  Improve  the  accuracy  of  the  calculation  of  Ez 
along  that  surface.  This  field  value  is  used  by  the  code  to  determine  the 
amount  of  charge  emitted  from  those  surface  cells  each  timestep.  The 
newly  emitted  electrons  must  furthermore  travel  some  finite  distance  from 
the  cathode  before  achieving  their  maximum  drift  energies.  Higher  than 
average  resolution  of  particle  motion  is  desired  over  this  critical  phase 
so  aZ  is  increased  only  gradually  over  the  next  1.43  cm.  Specifically, 
the  11  cells  from  Z  *  1.08  to  Z  =  2.51  cm  are  given  the  following  aZ 
values:  0.08,  0.09,  0.10,  0.11,  0.12,  0.13,  0.14,  0.15,  0.16,  0.17,  and 
0.18  cm.  Filling  the  bulk  of  the  tube  length,  from  Z  *  2.51  to  Z  *  23.79 
cm,  are  112  axial  cells  of  aZ  =  0.19  cm  each.  Thereafter,  covering  the 
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1.43  cm  to  the  tip  of  the  rod  cathode,  follows  another  set  of  11  cells  of 
Identical  values  but  reverse  order  compared  to  those  at  the  hollow  cathode 
tip.  From  that  point,  Z  *  25.22  cm,  up  to  the  right-hand  boundary  of  the 
computational  region  at  Z  *  28.26  cm.  Is  an  extension  of  the  vacuun  feed 
coaxial  power  transmission  line  to  the  rod  diode.  Since  no  particles  will 
flow  in  that  region,  spatial  resolution  may  be  relaxed.  Thus,  over  the 
next  1.12  cm,  the  axial  grid  size,  aZ,  Increases  in  the  order;  0.10,  0.12, 
0.14,  0.16,  0.18,  0.20,  and  0.22  cm;  In  7  cells.  Finally,  8  cells  of  a 
constant  aZ  «  0.24  cm  fill  the  remainder  of  the  computational  mesh.  Thus, 
there  are  a  total  of  157  cells  in  the  axial  dimension.  Adding  a  necessary 
monolayer  of  boundary  cells  in  each  dimension  yields  a  net  (158  x  25  =) 
3950  grid  points  in  the  entire  numerical  model.  To  complete  the  numerical 
model  of  the  CPA  on  this  grid,  the  graphite  hollow  beam  dump  is 
represented  by  a  perfectly  absorbing,  perfectly  conducting  solid,  filling 
the  cross-section  from  Z  =  25.22  to  Z  =  28.26  era  and  from  R  =  1.60  to  R  = 
2.40  cm.  It  is  considered  to  be  electrically  part  of  the  cylindrical 
anode  tube.  Finally,  the  central  rod  cathode  is  approximated  by  a  perfect 
conductor  of  radius,  0.32  cm,  lying  along  the  axis  from  Z  *  25.22  to  Z  * 
28.26  cm.  The  model  device  which  emerges  is  drawn  to  scale  in  Figure  3. 

Having  established  the  number  of  computational  cells,  an  estimate  can 
be  made  of  the  number  of  particles  needed  In  the  typical  simulation.  If 
the  two  beams  are  roughly  confined  to  only  the  high-resolution  radial 
regions,  they  will  normally  be  restricted  to  a  maximun  of  (6  +  8)  x  143  * 
2002  cells.  Requiring,  for  numerical  noise  prevention,  about  10  to  20 
macroelectrons  per  cell  yields  a  goal  of  20,000  to  40,000  particles  in 
each  simulation  where  both  beams  are  present.  These  electrons  are  emitted 
onto  the  grid  from  the  radial  faces  of  the  two  cathodes  according  to 
HAGIC's  Gaussian  emission  algorithm.  At  the  onset  of  a  run,  the  mesh  is 
devoid  of  particles.  Immediately,  the  hollow  cathode  is  "turned  on" 
through  the  imposition  of  a  short-rise  time  potential  pulse  of  -1.5  MV 
compared  with  the  0.0  volts  of  the  "grounded"  anode  drift  tube.  The  shape 
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of  this  pulse  Is  shown  In  Figure  4a.  The  hollow  beam  requires  at  least 
0.805  nanosecond  to  cross  to  the  "graphite"  beam  stop  at  the  other  end  of 
the  tube.  To  give  the  hollow  beam  time  to  equl liberate,  2.12  ns  (just 
over  1500  time  steps)  are  allowed  to  pass  before  a  voltage  pulse  Is 
applied  to  the  rod  cathode.  In  the  NRL  experiments,  the  voltage  was 
ramped  from  0  to  -1.0  MV  on  the  rod  cathode  over  about  5  nanoseconds.  For 
computational  efficiency,  this  pulse  was  ramped  in  two  linear  stages  to 
-1.0  MV  over  1.5  nanoseconds.  This  second  pulse  is  depicted  in  Figure 
4b.  There  is  first  a  rapid  (0.0925  ns)  rise  to  -0.5MV,  followed  by  a  more 
gradual  escalation  to  the  full  -1.0MV  potential.  Thereafter  another  1.0 
nanosecond  Is  necessary  for  the  1.0  MeV  electrons  to  cross  the  tube.  This 
yields  a  total  run  time  of  about  4.5  nanoseconds  or  3215  timesteps  of  At  s 
1.40  x  10“ 12  second  each  for  a  complete  simulation  with  both  beams 
present. 


» 


» 
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(b)  Rod  Cathode  Voltage. 


Figure  4. 


Voltage  Pulses  Applied  in  the  CPA  Simulations 
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4.  THE  MAGIC  SIMULATION  CODE 


MAGIC  (MAGnetic  2nsu^a^on  Code)  is  a  computer  code  designed  to 
simulate  self-consistent  (nonlinear)  processes  typically  associated  with 
pulsed  power  systems.  Examples  include  power  transport,  beam  propagation, 
magnetic  insulation,  diode  operation,  microwave  tube  design,  etc.  MAGIC 
is  an  intermediate  size  code  (about  20,000  statements),  and  is  highly 
optimized  both  for  user  efficiency  and  computational  efficiency.  It  was 
developed  by  MRC  under  contract  to  Sandia  National  Laboratories  for  use  in 
their  light  ion  and  electron  fusion  program.  MAGIC  is  now  entering  the 
fifth  year  of  development  as  a  production  code.  There  is  an  active  user 
group,  supported  by  an  extensive  User's  Manual  and  active  consultation 
with  the  authors.  The  code  is  available  on  Cyber  176,  CRAY-1,  and  VAX-780 
computers.  The  code  has  been  docunented  extensively  (References  4-6)  and 
many  results  have  been  published  (e.g..  References  7-9);  the  discussion 
which  follows  will  deal  with  the  code's  options  and  features  which  are 
crucial  to  successful  simulation  of  relativistic  electron  beam  devices. 

Conducting  surfaces  can  be  imposed  arbitrarily  to  exactly  duplicate 
structural  geometry  in  any  of  four  coordinate  systems  (cartesian, 
cylindrical  R-Z,  cylindrical  R-e,  and  spherical).  Periodic  and  mirror 
symmetry  conditions,  which  apply  to  both  fields  and  particles,  can  be 
imposed  on  the  ends  of  a  computational  region.  MAGIC  has  the  capability 
of  linking  one-dimensional  transmission  lines  to  the  two-dimensional 
simulation;  this  provides  the  means  to  model  the  interior  structure  of 
devices,  which  are  plasma  free.  Such  lines  can  be  truncated  with  a 
"lookback"  model  which  allows  reflected  waves  to  escape  from  the  system. 

Externally  applied  (static)  magnetic  fields  can  be  imposed  in  any 
direction  and  can  be  a  function  of  the  spatial  coordinates.  (The  dynamic 
magnetic  field  adds  linearly  to  the  static  field  to  affect  particle 
motion.)  The  functional  form  of  electric  fields  from  a  static  voltage  can 
be  determined  initially  from  the  solution  of  Poisson's  equation.  These 
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fields  can  be  imposed  by  means  of  an  external  circuit  to  represent  a  power 
pulse  propagating  in  the  third  dimension.  Either  the  steady-state  or  the 
time  dependent  case  can  be  represented.  In  either  case,  the  dynamic 
fields  again  add  linearly  to  the  static  fields  to  affect  particle  motion. 

The  emission  of  electrons  or  ions  can  be  simulated  over  arbitrary 
predetermined  conductor  surfaces  by  application  of  Poisson's  equation  to 
surface  cells,  i.e.,  cells  immediately  adjacent  to  the  electrode  surface. 
The  external  normal  electric  field  and  charge  existing  within  the  cell  are 
used  to  determine  the  surface  normal  field.  To  obey  Child-Langmuir 
physics,  this  field  should  vanish.  In  MAGIC,  enough  charge  is  created 
near  the  surface  to  force  this  field  to  vanish  (if  the  field  has  the  wrong 
sign,  no  charge  is  created).  This  charge  is  normally  broken  up  into  a 
small  number  of  macroparticles,  which  are  then  distributed  randomly  along 
the  cell  surface.  A  small  spatial  distribution  in  the  normal  direction 
provides  effective  temporal  separation.  The  surface  force  algorithm  is 
also  altered  to  take  actual  surface  charge  into  account. 

This  very  simple  model  has  seen  effective  use  in  a  variety  of 
applications.  It  has  been  validated  by  comparison  with  one-diiw>«< ional 
Child-Langmuir  analysis.  The  primary  requirement  to  achieve  accuracy  is 
fine  normal  spacing  near  the  surface,  and  this  is  easily  attained  in  MAGIC 
using  the  nonuniform  spacing  option.  Finally,  MAGIC  has  the  capability  to 
simultaneously  model  up  to  10  particle  species,  which  may  be  useful  if 
positive  ions  are  believed  to  play  any  role. 

Two  standard  explicit  (centered-difference  and  time- symmetric)  and  one 
novel  implicit  (time-biased)  electromagnetic  field  solving  algorithms  are 
available  as  user  options  in  the  code.  Particular  mention  should  be  made 
of  the  implicit  solver  since  this  is  the  one  actually  employed  in  these 
simulations.  The  basic  idea  behind  this  algorithm  is  very  similar  to  the 
notion  behind  predictor-corrector.  That  is,  spatial  gradients  of  a 
magnetic  field  at  several  points  in  time  are  used  to  estimate  the  time 


rate  of  change  of  the  electric  field.  If  these  several  points  in  time 
include  advanced  (unknown)  magnetic  fields,  then  the  algorithm  is 
implicit.  In  MAGIC,  the  implicit  equations  are  solved  by  repetition  for  a 
single  time  step,  making  use  of  relaxation  coefficients  to  improve 
convergence.  This  algorithm  was  originally  developed  by  B.  Godfrey  of  MRC 
(Reference  10),  and  has  been  extensively  tested  by  use  (e.g. ,  Reference  6) 
in  applications  involving  magnetic  insulation. 

The  primary  advantages  of  the  time-biased  over  standard  centered 
difference  methods  are  that  (1)  high  frequency  noise  is  greatly  reduced, 
(2)  certain  numerical  instabilities  are  damped,  and  (3)  the  Courant 
criterion  is  usually  relaxed.  On  this  last  point,  the  penalty  of 
iterative  solution  outweighs  the  advantage  of  larger  time  steps.  That  is, 
the  time-biased  algorithm  is  always  more  expensive  computationally. 
However,  use  of  this  algorithm  (along  with  temporal  filtering  of  magnetic 
forces)  has  virtually  eliminated  the  debilitating  wave  trapping  normally 
encountered  in  magnetic  insulation  type  problems  (see  Reference  7  for  a 
discussion  of  wave  trapping). 

Local  charge  conservation  is  another  algorithm  which  is  essential  to  a 
good  representation  of  intense  particle  beam  devices.  The  majority  of  PIC 
codes  conserve  charge  only  in  the  global  sense  -  that  is,  system  charge 
equals  charge  created  minus  charge  destroyed.  In  many  applications,  the 
fluctuations  in  individual  cell  fields  (due  to  various  charge  weighting 
schemes)  are  small  compared  with  fields  themselves.  In  such  applications, 
the  global  conservation  schemes  are  adequate.  In  beam-diode  problems, 
however,  the  accumulated  charge  creation  and  destruction  typically  greatly 
exceed  the  system  charge.  In  this  case,  the  field  fluctuations  can  exceed 
the  magnitude  of  the  fields  themselves.  This  proceeds  linearly  with  time 
in  a  simulation,  culminating  in  obviously  catastrophic  results. 

This  effect  is  eliminated  if  charge  conservation  is  ensured  locally 
(cell -by-cel  1 ) .  There  are  several  approaches  possible  to  achieve  this. 


For  example,  some  electromagnetic  PIC  codes  periodically  during  the 
simulation  (say,  every  50  time  steps)  will  solve  Poisson's  equation,  and 
use  the  results  to  "correct"  the  local  conservation  directly  through  the 
current  density  algorithm.  By  contrast,  MAGIC  relies  on  consideration  of 
multiple,  orthogonal  motions  of  particles  and  exact  enforcement  of  the 
continuity  equation  in  the  direction  of  motion.  Thus,  local  conservation 
is  provided  automatically  by  the  current  density  algorithm  in  MAGIC, 
eliminating  the  need  for  subsequent,  periodic  "corrections". 

MAGIC  has  been  used  to  simulate  too  many  devices  to  conveniently  list 
them  all  here.  Still,  it  is  useful  to  mention  several  examples  which 
illustrate  the  power  and  feasibility  of  the  code.  The  first  is  explained 
in  Reference  11.  MAGIC  was  there  used  to  simulate  the  effect  on  electron 
beam  transmission  of  adding  a  dielectric  liner  to  the  cylindrical  outer 
wall  of  an  experimental  drift  tube  cavity.  The  second  example  relates  to 
an  extensive  series  of  foil less  diode  simulations  conducted  with  MAGIC  by 
Sandia  National  Laboratories  (Reference  12).  Some  22,000  macroelectrons 
were  used  in  these  runs  which  tested  the  effects  of  different  values  of  an 
applied,  uniform,  axial  magnetic  field,  B0.  The  final  example  presented 
in  Reference  13  looked  at  the  problem  of  electron  leakage  in  the  NRL 
hybrid-inverted  coaxial  magnetron. 


5.  RESULTS 


The  primary  results  from  this  research  effort  were  derived  from  five 
separate  sets  of  MAGIC  simulations  conducted  on  the  NRL  VAX-11/780 
computer.  As  detailed  in  the  following  five  subsections,  each  simulation 
shed  light  on  different  aspects  of  the  experimental  device  being  modeled. 
They  represent  a  logical  progression  of  added  insight.  As  new  data 
emerged  from  each  run,  minor  changes  were  made  to  the  original,  numerical 
model  and  new  diagnostic  output  was  requested.  The  starting  point  of  the 
effort  was  a  simulation  of  the  full  device  as  shown  in  Figure  2  with 
voltage  pulses  applied  as  shown  in  Figure  4  but  with  no  particle 
emission.  Then  the  simulation  was  rerun  with  only  the  hollow  cathode 
pulse  applied,  allowing  particle  emission  there,  to  test  the  formation  of 
the  outer  hollow  beam.  In  the  next  case,  both  pulses  were  applied  again 
and  both  cathodes  were  allowed  to  emit.  This  was  the  first  test  of  the 
dual  beam  configuration  and  it  appeared  to  confirm  previous  experimental 
results.  To  gain  additional  physical  information,  the  central  beam  was 
then  tested  without  the  hollow  beam  present.  Finally,  the  dual  beam  case 
was  rerun,  this  time  using  a  shortened  voltage  pulse  time  history  for  the 
rod  cathode.  The  major  conclusions  to  be  drawn  from  all  of  these  runs  are 
summarized  in  Section  6. 

5.1  "Cold  Test"  of  the  CPA.  A  useful  first  test  of  the  CPA  configuration 
being  modeled  was  a  simple  "cold  test"  of  the  device.  This  involved 
imposing  the  desired  voltage  pulses  (see  Figure  4)  on  the  respective, 
opposing  diode  ports  while  not  allowing  any  particle  emission  from  the  two 
cathodes.  This  showed  how  the  two  initially  pure  TEM  waveforms  combine  in 
the  cavity  and  Illustrated  the  nature  of  the  net  electric  fields  down  the 
entire  length  of  the  source- free  tube. 

At  the  start  of  the  simulation,  a  voltage  of  1.5  x  106  volts  was 
applied  with  a  rise-time  of  less  than  0.1  ns  across  the  entrance  port  to 
the  hollow  cathode  from  its  power  transmission  line.  This  was  introduced 
as  a  pure  TEM  wave  at  that  numerical  boundary  with  a  radially  graded 


electric  field  of  about  (1.5  x  106)/(0.38  x  1Q~2)  *  -3.95  x  108  volts  per 
meter.  At  this  point  in  the  presentation,  a  certain  feature  of  the  MAGIC 
code's  treatment  of  boundary  conditions  should  be  explained.  The  same 
input  data  prescription  that  allows  the  imposition  of  the  TEM  wave  along 
the  specified  numerical  boundary  segment  also  ensures  the  absorption  of 
any  TEM  waves  incident  on  that  surface.  This  eliminates  unwanted 
reflections  of  waves  from  that  boundary  as  long  as  the  incident  wave  is 
purely  TEM.  To  guarantee  that  as  closely  as  possible,  the  axial  lengths 
of  each  of  the  two  transmission  line  entrance  port  regions  in  the 
numerical  model  should  be  much  greater  than  the  respective  radial 
separations  between  the  coaxial  conductors  in  those  lines.  For  the  hollow 
cathode  port,  the  axial  length  is  1.08  cm  (see  Figure  2)  compared  to  the 
radial  gap  of  0.38  cm,  thus  yielding  a  comfortable  ratio  of  better  than 
2.8:1.  For  the  rod  cathode  inlet,  the  corresponding  lengths  for  this 
first  simulation  were  2.56  cm  compared  to  1.28  cm  or  a  ratio  of  exactly 
2:1  which  was  not  quite  as  comfortable. 

After  the  passage  of  500  time  steps  of  1.4  x  10" 12  second  each,  the 
wave  front  had  traveled  about  21  cm  down  the  tube  and,  thus,  had  almost 
reached  the  opposing  diode  structure.  It  was.  decided  to  sample  the 
electric  field  structure  in  the  tube  at  that  time.  Figure  5  displays  a 
MAGIC  "LOGPRINT"  map  of  the  radial  electric  field  component  value  at  every 
data  point  in  the  mesh.  A  word  of  explanation  regarding  the  plot  is  in 
order.  The  upper  half  of  the  figure  represents  that  half  of  the  device 
closer  to  the  rod  diode  and  the  lower  half  represents  the  half  closer  to 
the  hollow  cathode,  respectively.  In  each  half,  the  lower  margin  consists 
of  a  column  of  1-  to  3-digit  numbers  above  which  is  a  blank  column.  That 
blank  column  corresponds  to  the  tube  centerline.  The  1-  to  3-  digit 
numbers  immediately  below  it  are  the  axial  data  cell  position  numbers  (see 
Figure  2).  The  left  margin  of  the  upper  half  of  the  figure  consists  of 
two  vertical  columns.  The  first  is  the  "tens"  digit  and  the  second  is  the 
"units"  digit  of  the  radi al  data  cell  position  numbers.  In  the  overall 
map,  the  silhouettes  of  the  two  cathodes  and  the  anode  beam  stop  are 
clearly  visable.  Note,  however,  that  their  respective  left-to-right 
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Figure  5.  Er  at  t  =  500 At  for  "Cold  Test"  Case 


positions  have  been  reversed  compared  to  the  cross  section  portrayed  in 
Figure  2.  The  numbers  which  appear  at  each  data  point  in  Figure  5 
represent  the  order-of-magnitude  (power  of  ten)  of  the  respective  value  of 
Er  (in  volts  per  meter)  with  a  scale  factor  of  1.  Minus  signs  are 
overprinted  on  any  negative  value  of  Er  on  the  mesh.  Stand  alone  minus 
signs  indicate  unprintably  small  negative  values.  The  applied  radial 
electric  field  values  of  order  10 8  V/m  are  clearly  present  at  the  hollow 
cathode  power  input  port.  The  field  strengths  drop  off,  as  expected,  when 
one  scans  leftward  toward  the  wave  front.  The  traveling  front  is  clearly 
outlined  at  around  z  =  132.  The  corresponding  "LOGPRINT"  map  for  the 
axial  component  of  the  electric  field  at  the  same  instant  in  time  is  shown 
in  Figure  6.  Note  that  Ez  is  everywhere  negative  as  expected. 

Since  it  was  not  necessary  in  this  "cold  test"  to  worry  about  a  hollow 
electron  beam  being  formed  and  then  reaching  equilibrium,  the  turn-on  of 
the  rod  cathode  voltage  pulse  was  advanced  from  the  t  =  2.12  ns  shown  in 
Figure  4  to  t  =  1.5  ns  instead.  The  shape  of  the  pulse  itself,  however, 
was  kept  the  same.  By  t  =  3.5  ns  (i.e.,  t  «  2500  at),  an  approximate 
equilibrium  of  the  fields  in  the  device  should  have  been  approached. 
Figures  7  and  8  show  the  “LOGPRINTs"  of  the  radial  and  axial  electric 
field  components,  respectively  at  that  time.  Some  minor,  low  magnitude 
ripples  are  still  apparent  in  the  Er  map  but  that  for  Ez  is  much 
smoother.  The  axial  field  is  predominantly  negative  from  the  hollow 
cathode  all  the  way  out  to  about  z  =  90  with  the  remainder  turned  positive 
due  to  the  opposing  rod  cathode  voltage  pulse.  The  field  strengths  are 
particularly  high  near  the  tips  of  the  cathodes  as  expected  due  to  two- 
dimensional  field  enhancement  effects  (Reference  14).  It  can  further  be 
seen  that  Ez  along  the  hollow  cathode  power  input  boundary  is  on  the 
order  of  104  V/m  compared  to  Er  values  of  108  there.  This  implies  an 
almost  perfect  TEM  field  structure  there  which,  in  turn,  guarantees  almost 
no  false  reflections  there.  The  situation  is  not  quite  so  favorable  along 
the  rod  cathode  power  inlet  at  the  opposite  end  of  the  tube.  There,  Ez 
is  on  the  order  of  1Q5-106  V/m  while  Er  is  like  10M08  V/m.  Thus,  less 
than  two  orders  of  magnitude  separate  the  two  components  there. 
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Figure  6.  Ez  at  t  *  500  At  for  "Cold  Test"  Case 
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Figure  7.  Er  at  t  »  2500  &t  for  "Cold  Test"  Case 
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Figure  8.  Ez  at  t  =  2500  At  for  "Cold  Test"  Case. 


Implying  deviations  from  a  pure  TEM  waveform  of  at  least  several  percent. 
Some  low-order  but  noticeable  reflections  should  take  place  there.  More 
on  that  later. 

The  next  important  diagnostic  output  from  this  run  consisted  of  axial 
profiles  of  the  electric  field  components  along  the  paths  which  would  be 
traveled  by  the  two  electron  beams.  These  profiles  are  plotted  in  Figure 
9.  The  upper  two  profiles,  taken  in  combination,  plot  Ez  (in  V/m) 
versus  z  (in  meters)  across  the  full  length  of  the  device  for  a  constant 
tube  radius  of  r  =  1.94  cm,  which  should  be  about  0.08  cm  less  than  the 
average  outer  radius  of  the  hollow  electron  beam.  Note  that  the  field 
strength  scales  of  the  two  plots  are  different;  the  left  one  runs  from  -8 
x  108  to  2  x  108  V/m  while  the  right  has  values  only  from  -1.5  x  107  to 
0.5  x  107  V/m.  The  plots  show  a  strong,  negative  Ez  over  the  first  2 
cm  out  from  the  hollow  cathode  tip  with  a  complimentary  2  cm  length  of 
enhanced  Ez  at  the  opposite  end  of  the  tube,  just  before  the  anode  beam 
stop.  The  peak  magnitude  of  Ez  at  that  far  end  is  down  by  a  factor  of 
about  50  from  that  near  the  cathode.  These  plots  show  in  much  stronger 
terms  than  Figure  8,  that  Ez  is  effectively  negligible  over  most  of  the 
tube's  length  except  at  the  ends. 

The  bottom  two  plots  of  Figure  9  give  the  axial  profiles  of  Er 
(left)  and  Ez  (right)  over  most  of  the  tube's  length  at  a  constant 
radius  r  *  0.24  cm  which  is  about  0.08  cm  inside  the  expected  outer  radius 
of  the  central  electron  beam  envelope.  The  radial  field  component  plot  is 
quite  unremarkable.  The  strong  Er  near  the  tip  of  the  rod  cathode  can 
be  expected  to  impart  a  significant  radial  momentum  component  to  the 
motion  of  electrons  emitted  there  but  this  is  to  be  expected  for  both 
beams.  More  significant  is  the  profile  for  Ez  at  this  radius.  Near  the 
rod  cathode  tip,  a  strong,  positive  Ez  will  accelerate  the  electrons 
down  the  tube  centerline  but  at  the  opposite  end  of  the  tube,  a  region  of 
strong  Ez  of  opposite  sign  is  encountered.  There,  the  central  beam 
electrons  will  be  decelerated.  This  is  to  be  expected  since  the  hollow 
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Figure  9.  Axial  profiles  of  Ez  at  R  =.20  (top  two) 
Ez  at  R  =  5  (bottom  right)  and  Er  at 
R  =  5  (bottom  left). 


cathode  is  at  -1.5MV  while  the  rod  cathode  is  at  only  -1.0MV.  In  the 
actual  experiment,  the  tube  length  of  over  5  meters  makes  electron 
reflection  from  such  a  reversed  field  region  ignorable  when  considering 
early  time  (<10  ns)  results  since  the  central  electrons  will  require  at 
least  15  ns  just  to  travel  once  to  the  region  of  the  hollow  cathode.  In 
these  simulations,  however,  the  shorter  tube  length  of  only  25  cm  requires 
that  care  be  taken  to  sort  out  this  reflecting  electron  effect  from  the 
gross  behavior  of  the  central  beam.  That  question  is  addressed  further  in 
subsections  5.3,  5.4,  and  5.5. 

The  final  major  diagnostic  output  from  the  code  run  consisted  of  the 
three  time  history  plots  shown  in  Figure  10.  Examining  them  from  left  to 
right,  the  first  presents  the  temporal  shape  of  the  actual  voltage 
measured  radially  between  the  hollow  cathode  and  the  anode  tube.  The  TEM 
wave  applied  at  the  power  pulse  inlet  boundary  was  only  1.5MV  (peak).  The 
non-emitting  hollow  cathode  tip,  however,  represents  an  infinite  impedance 
mismatch  in  the  coaxial  feed  line  and  therefore  reflects  the  input  wave 
and  doubles  the  resultant  hollow  cathode-to-anode  peak  voltage  to  3.0MV. 
It  is  this  voltage  that  would  be  measured  experimentally.  If  the  hollow 
cathode  emits  electrons,  that  changes  the  impedance  mismatch  which  changes 
the  wave  reflection  coefficient  which  changes  the  measured  voltage.  With 
the  hollow  beam  present,  this  was  experimentally  measured  to  be  1.5MV. 
Therefore,  the  voltage  applied  to  that  boundary  in  the  simulation  must  be 
adjusted  downwards  until  it  yields  a  net  measured  electrode  voltage  of 
1.5MV.  This  was,  in  fact,  done  as  additional  simulations  were 

accompl i shed. 

The  next  plot  in  Figure  10  is  the  corresponding  measured  voltage  pulse 
on  the  rod  cathode  with  respect  to  the  anode  (specifically,  the  anode  beam 
dump).  Comparing  this  and  the  previous  one  to  their  ideal,  applied, 
numerical  counterparts  as  shown  in  Figure  4,  demonstrates  a  very  close 
correspondence.  The  hollow  cathode  pulse  is  nearly  identical.  The  rod 
cathode  pulse  (ignoring  the  previously  mentioned  turn-on  time  advance) 
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also  has  good  shape  fidelity  except  for  the  2-5%  ripples  which  appear  to 
be  due  to  minor  reflections  of  waves  hitting  the  rod  cathode  power  inlet 
boundary.  (It  was  already  stated  that  the  non-TEM  components  will  be 
reflected  there.)  To  help  reduce  these  ripples,  the  axial  length  of  the 
coaxial  power  inlet  in  the  model  was  subsequently  lengthened  by  the 
addition  of  two  axial  cell  columns  of  tH  =  0.24  cm.  That  addition  brings 
the  total  axial  cell  count  up  from  156  to  158  cells,  giving  a  total  of 
3950  data  cells  covering  the  entire  numerical  mesh. 

The  final  plot  in  Figure  10  simply  presents  the  time  history  of  the 
magnitude  of  Ez  (in  volts/meter)  measured  just  opposite  the  tip  of  the 
rod  cathode.  It  faithfully  reproduces  the  temporal  shape  of  the  rod 
cathode  voltage  pulse.  No  interference  from  the  opposing  cathode  voltage 
pulse  is  obvious  from  this  plot.  Without  charged  particles  present  to 
strongly  couple  them,  the  two  pulses  from  either  end  of  the  accelerator 
tube  seem  fairly  independent.  No  threat  to  the  successful  operation  of 
this  device  can  yet  be  observed.  The  next  step  in  these  studies  involved 
the  "turning  on"  of  emission  from  the  hollow  cathode  tip  alone  and 
examining  the  formation  of  the  hollow  electron  beam.  In  the  CPA 
experiments  at  NRL,  no  difficulties  were  ever  experienced  with  propagating 
a  hollow  beam  of  some  30  kA  down  the  full  length  of  the  five  meter  tube 
whether  or  not  the  central  beam  was  present. 

5.2  Hollow  Beam  Only.  Haying  tested  the  pure  electromagnetic  response  of  , 
the  CPA  model,  the  next  step  involved  the  inclusion  of  charged  particles 
in  the  simulations.  In  the  CPA  experiments,  it  was  observed  that  the 
hollow,  high  current,  electron  b«am  would  always  faithfully  propagate  down 
the  entire  length  of  the  tube.  In  this  first  numerical  run  with 
particles,  it  was  decided  to  examine  the  formation  of  the  hollow  beam 
alone,  without  any  voltage  pulse  being  applied  to  the  rod  cathode  and 
without  electron  emission  being  allowed  there.  For  this  simulation,  a 
peak  voltage,  1.5MV  pulse  was  applied  to  the  hollow  cathode  diode.  The 
pulse  shape  was  that  shown  in  Figure  4a.  with  a  rise  time  of  slightly  less 


than  0.1  ns.  Electron  emission  was  permitted  along  the  entire  protruding 
front  face  of  the  hollow  cathode.  (Note  that  MAGIC  uses  a  physically 

•  true,  Gaussian  emission  algorithm  which  calculates  how  much  charge  should 
be  emitted  from  a  given  surface  at  a  given  time  step  by  sampling  the 
local,  normal  electric  field  -  see  Section  4  of  this  report.)  If  the 
hollow  cathode  diode  being  modeled  corresponded  exactly  to  the  physical 

•  device  under  study,  it  would  be  mandatory  to  include  in  the  simulation 
electrons  emitted  from  the  cathode  shank.  Such  shank  emission  has  been 
shown  to  significantly  effect  the  current  characteristics  of  a  diode 
(Reference  15).  Such  emission  was  not  included  here  for  the  following 

•  reasons: 


a.  The  radial  anode-cathode  gap  in  the  model  is  artificially  small 
and  would  therefore  give  rise  to  an  artificially  large  shank 
emission  current. 

b.  As  long  as  the  hollow  beam  is  stable  and  remains  within  reason¬ 
able  limits  with  regard  to  its  total  current,  any  finer  details 
with  regard  to  its  microscopic  characteristics  are  not  of  major 
importance  to  these  simulations. 

c.  The  "turning  on"  of  some  length  of  the  cathode  shank  to  emission 
would  require  a  lengthening  of  that  coaxial,  inlet  port  region 
in  he  numerical  model  in  order  to  ensure  accurate  TEM  wave  power 
injection  at  that  boundary.  This  would  require  the  addition  of 
extra  cells  to  the  numerical  model  which  would  increase  computer 
run  times.  (The  additional,  numerical  particles  emitted  from 
the  shank  would  even  more  dramatically  increase  typical,  run 
times.) 


After  500  time  steps  of  At  =  1.4  x  10" 12  second,  a  sample  electron 
position  plot  was  generated,  as  shown  in  Figure  11.  Strong  ripples  are 
evident  in  the  beam  as  is  typical  in  foil  less  diode  operation  (Reference 
16).  It  appears  that  substantial  phase  mixing  of  the  ripples  takes  place 
by  halfway  down  the  tube.  In  viewing  this  figure,  it  is  important  to  keep 
in  mind  the  gross  distortion  of  axial -to-radial  scales.  The  total 
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distance  from  the  dashed  centerline  to  the  outer  tube  wall  is  actually 
only  2.4  cm  while  the  distance  from  the  hollow  cathode  tip  to  the  axially- 
opposing  anode  beam  stop  is  24.14  cm.  Thus  axial  dimensions  are  extremely 
compressed.  This  effect  accentuates  the  amplitude  of  the  ripples. 

After  another  500  time  steps,  the  beam  is  fairly  well  established  as 
can  be  seen  from  the  sample  position  plot  of  Figure  12.  The  phase  mixing 
of  the  ripples  is  even  more  pronounced  here.  Significant  gaps  between  the 
particle  streams  near  the  cathode,  however,  indicate  the  advisability  of 
increasing  the  total  number  of  particles  used  to  represent  the  beam.  An 
increased  number  density  should  also  result  in  decreased  choppiness  of  the 
downstream  portion  of  the  beam.  In  this  simulation,  the  total,  average 
particle  count  was  about  19,000.  The  field  characteristics  of  this  single 
beam  model  were  likewise  diagnosed  in  detail. 

The  “LOGPRINTs"  for  both  the  radial  and  axial  components  of  the 
electric  field  are  displayed  in  Figures  13  and  14  respectively.  The  Er 
plot  shows  the  expected  maintenance  of  a  strong,  uniform  radial  E-field 
between  the  beam  and  the  anode  wall.  The  corresponding  E2  plot  shows 
high  field  strength  values  concentrated  at  either  end  of  the  tube  with  low 
values  in  between.  To  help  visualize  what  has  happened  to  the  E-field  in 
the  tube,  it  would  be  instructive  to  view  an  equipotential  plot.  However, 
MAGIC  works  only  with  Er( z,r)  and  Ez(z,r),  not  with  ^z.r).  Never¬ 
theless,  parallel  results  are  available  from  another  foil  less  diode 
simulation  the  author  has  run  (Reference  16)  and  these  are  shown  in  Figure 
15.  This  simulation  was  run  with  the  quasistatic  DI0DE2D  code  at  the 
Naval  Research  Laboratory.  The  top  part  of  the  figure  shows  the  numerical 
model  and  the  bottom  shows  <j»(z,r)  in  the  anode  tube  when  a  well -colli mated 
hollow  electron  beam  is  propagating  down  the  tube.  Note  that,  as  in  the 
present  simulation,  electrons  emitted  from  the  cathode  tip  are  first 
accelerated  axially  in  a  region  of  enhanced  Ez.  That  first  acceleration 
gradient  does  not  impart  the  full  diode  voltage  to  them.  For  the  example 
shown,  ony  four  10%  equipotential  contours  are  crossed  in  that  first 
region  before  electrons  effectively  enter  a  "drift"  mode  down  most  of  the 
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Figure  13.  Er  at  t  »  1500 at  for  Hollow  Beam  Case 
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remaining  length  of  the  tube.  Finally,  just  before  impacting  the  anode 
end-plate,  the  electrons  cross  another  region  of  enhanced  Ez  over  which 
they  are  brought  up  to  full  diode  energy.  Notice  that  the  regions  of 
enhanced  Ez  at  both  ends  of  the  beam  extend  radially  inward  with  only 
slightly  decreasing  intensity.  Carrying  these  qualitative  results  back  to 
the  present  MAGIC  simulation,  implies  that  the  presence  of  the  intense, 
outer,  hollow  beam  would  create  a  region  of  enhanced  Ez  of  the  wrong 
sign  at  the  tip  of  the  rod  cathode  which  would  have  to  be  overcome  by  the 
1.0MV  voltage  pulse  there  for  axial  propagation  of  the  central  beam.  This 
may  be  surmised  as  the  major  cause  for  the  propagation  problems 
encountered  in  the  actual  experiments  and  the  remainder  of  these  simula¬ 
tions  give  confirmation  to  that  hypothesis. 

In  agreement  with  the  above  description  for  a  typical  electron 
acceleration  history  in  the  hollow  beam  are  the  axial  profiles  of  Ez  at 
r  =  1.94  cm  (just  inside  the  outer  radius  of  the  hollow  beam)  displayed  in 
the  top  two  plots  of  Figure  16.  Large  values  of  Ez  are  axially  confined 
to  narrow  regions  (<2  cm)  immediately  adjacent  to  the  cathode  tip  (left) 
and  the  anode  beam  stop  (right).  The  bottom  two  plots  are  the 
corresponding  axial  profile  for  the  radial  electric  field  at  r  =  1.94  cm. 
They  show  nothing  particularly  remarkable  and  are  included  only  for  the 
sake  of  completeness.  Much  more  significant  is  the  axial  profile  of  Ez 
at  r  *  0.24  cm  (just  inside  the  future  outer  envelope  of  the  central  beam) 
shown  in  Figure  17.  As  could  be  predicted,  a  strong,  negative  Ez 
dominates  the  region  near  the  rod  cathode  tip  (right  plot).  The  peak 
field  there  is  -1.2MV/cm.  That  must  be  neutralized  by  the  rod  cathode 
voltage  before  electrons  can  even  be  emitted,  let  alone  accelerated  down 
the  axis.  In  the  opposite  half  of  the  device  (left  plot)  at  the  same 
radius.  Figure  17  shows  two  rather  extended  regions  of  significant 
negative  Ez  which  would  act  to  decelerate  any  central  beam  electron. 
The  simulations  in  the  following  three  sections  will  demonstrate  how  this 
field  structure  is  modified  by  the  presence  of  the  central  beam  and  how 
that  beam  is  effected  by  the  resultant  fields. 


Figure  17.  Axial  Profile  of  Ez  at  R  =  0*24  cm  and 
t  =  1000  At  for  Hollow  Beam  Case. 


The  final  set  of  diagnostics  for  this  Hollow-Beam-Only  case  is 
presented  in  Figure  18.  The  upper  left  plot  shows  the  measured  time 
history  <*f  the  voltage  between  the  hollow  cathode  and  the  anode.  Just  as 
in  the  previous  "cold  test"  simulation,  a  voltage  of  1.5  MV  is  being 
applied  at  the  inlet  boundary  but  is  then  being  augmented  by  the  wave 
reflected  from  the  impedance  mismatch  at  the  end  of  the  cathode  shank. 
With  no  hollow  beam,  the  mismatch  was  virtually  infinite,  resulting  in  a 
voltage  doubling.  Now,  with  the  hollow  beam  present,  there  is  only  a 
finite  mismatch  which  boosts  the  measured  peak  voltage  from  the  applied 
1.5  MV  to  about  2.2  MV.  Therefore,  to  agree  with  the  experiment,  the 
applied  voltage  must  be  lowered  to  about  (1.5)  (1.51/2.2  *  1.0  MV.  That 
is  the  value  used  there  in  all  of  the  following  simulations.  The  upper 
right  plot  monitors  the  magnitude  of  the  axial  electric  field  component  a 
distance  of  0.04  cm  axially  away  from  the  rod  cathode  tip.  It  shows  the 
creation  of  a  strong  negative  field  there  due  to  the  presence  of  the 
hollow  electron  beam.  This  should  inhibit  the  emission  and/or  propagation 
of  an  electron  beam  from  that  cathode.  The  bottom  two  plots  of  Figure  18 
correspond  to  the  Be  probe- type  measurements  such  as  those  actually 
performed  in  the  experiment  (although  the  experimental  probes  displayed 
time- integrated  results).  The  left  plot  traced  the  value  of  B0  at  a 
radius  of  2.345  cm  and  a  distance  4.185  cm  axially  downstream  from  the 
hollow  cathode  tip.  It  shows  a  rapid  leveling  off  of  BQ  to  an  average 
value  of  about  -0.675  webers  per  meter2  (  *  -6.75  kilogauss).  This 
translates  to  a  beam  current  of  about  79  kA  which  is  more  than  a  factor  of 
two  above  that  used  in  the  experiments.  Lowering  the  effective  diode 
voltage  from  2.2  MV  to  1.5  MV  in  the  remaining  simulations  lowered  the 
hollow  beam  current  to  a  more  reasonable  value.  The  lower  right  plot  of 
Figure  18  traces  B0  at  the  same  radius  but  at  an  axial  distance  of 
20.395  cm  downstream  from  the  cathode  tip.  It  shows  the  beam 

time-of- flight  delay  and  the  leveling  off  to  the  same  value  of  about 
-0.675  w/m2. 
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Figure  18.  Hollow  Beam  Case  Time  Histories  for  (a) 
Hollow  Cathode  Voltage  Pulse,  (b)  Ez  at 
(142,5),  (c)  Bfi  at  (34,24).  and  (d)  Bfl 
at  (1^9,24)  0 


5.3  Dual  Beam  Case  I.  The  next  step  in  this  simulation  program  was  to 
test  the  dual  beam  operation  of  the  CPA  model  with  an  appropriately 
adjusted  hollow  cathode  voltage  value  and,  thereby,  a  lowered  hollow  beam 
current.  To  accomplish  this,  the  plateau  value  of  the  peak  voltage  in 
Figure  4a  was  lowered  to  1.0MV.  The  corresponding  peak  voltage  for  the 
rod  cathode  (Figure  4b)  was  left  as  shown.  Figures  19,  20,  and  21  display 
six  sample  particle  plots  for  the  electrons  of  both  beams  at  0.7,  2.1, 
2.8,  3.5,  4.2  and  4.9  nanoseconds.  (The  rod  cathode  voltage  pulse  began 
at  2.12  ns.)  The  peak  simulation  particle  count  for  this  run  was  about 
35,000  macroparticles.  The  increased  particle  density  significantly 
enhanced  the  smooth  appearance  of  the  equilibrium  hollow  beam.  For  the 
parameters  chosen,  these  figures  clearly  demonstrate  that  the  central  beam 
does  propagate  across  the  entire  CPA  tube.  That  is  qualitatively 
obvious.  To  quantify  the  nature  of  this  propagation,  more  detailed 
diagnostics  must  be  examined. 

The  fi.st  of  these  is  the  Ez-field  LOGPRINT  of  Figure  22.  Unlike 
the  corresponding  plots  of  Figures  5  and  6,  the  minus  signs,  which 
indicate  negative  field  values,  are  not  overprinted  on  the  field  magnitude 
number  but  rather  are  offset  printed  in  a  one-to-one  map  directly  above 
the  number  plot.  This  makes  the  negative  values  a  good  deal  easier  to 
pick  out.  As  in  the  previous  case,  the  highest  field  values  (7's  and  8's) 
are  most  strongly  concentrated  at  either  end  of  the  tube.  For  the  hollow 
beam,  there  are  intense  negative  Ez  regions  at  the  beginning  and  end 
while  those  for  the  central  beam  are  appropriately  positive  and  of  smaller 
extent.  Unlike  the  hollow  beam,  the  central  beam  also  hosts  regions  of 
7-strength  Ez  fields  located  far  from  its  end-points  (e.g.,  -  between  Z 
=  38  and  Z  =  57  in  Figure  22).  What's  more,  these  regions  are  of  mixed 
sign  which  indicates  some  degree  of  significant  beam  deceleration. 

The  corresponding  radial  E-field  LOGPRINT  may  be  found  in  Figure  23. 
As  could  be  expected,  the  field  is  almost  everywhere  negative  and  displays 
a  strong  uniformity  in  the  axial  dimension  due  to  the  axially  uniform 
electron  spau  charge  distribution.  Particularly  noteworthy  is  the  very 


Figure  19.  Hollow  Beam  at  t  =  500at  and  t  =  1500at  for  Dual  Beam  Case 
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Figure  20.  Both  Beams  at  t  =  2000< 


Figure  21.  Both  Beams  at  t  =  3000 at  and  t  =  3500 at  for  Dual  Beam  Case 
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Figure  23.  Er  at  t  =  3500  4t  for  Dual  Beam  Case  I. 
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uniform  region  of  8-strength  Er  sandwiched  between  the  hollow  beam  and 
the  anode  tube.  This  Is  Indicative  of  the  near  textbook  performance  of 
that  beam.  The  picture  of  Its  steady  behavior  is  further  enhanced  by  the 
axial  profiles  of  Ez  and  Er  as  measured  at  t  =  4.9  ns  at  a  radius  of 
r  *  1.94  cm.  (just  inside  the  hollow  beam's  outer  radius)  shown  in  Figure 
24.  Ez  in  the  beam  oscillates  mildly  about  zero  throughout  the  entire 
beam  except  in  the  end  regions  where  it  goes  strongly  negative. 
Conversely,  Er  maintains  a  rough  plateau  of  approximately  -1.10  x  108 
V/m  along  the  entire  beam  except  for  the  ends  where  it  goes  to  zero.  In 
contrast  to  these  rather  classic  axial  field  profiles  is  the  axial  profile 
of  Ez  at  a  radius  of  r  *  0.24  cm  (just  inside  the  outer  radius  of  the 
central  beam)  shown  at  the  same  time  in  the  simulation  in  Figure  25.  In 
the  half  of  the  tube  nearest  the  rod  cathode  (the  right-hand  figure)  the 
field  behaves  as  expected.  In  the  downstream  half  of  the  tube,  however,  a 
region  of  strong  deceleration  is  encountered  at  about  Z  =  0.09  m.  and  the 
final  peak  acceleration  field  is  more  than  a  factor  of  two  lower  than  that 
near  the  cathode. 

The  macroscopic  behavior  of  this  dual  beam  case  is  summarized  in 

Figure  26.  The  upper  two  plots  show  the  time  histories  of  the  hollow 

cathode  and  rod  cathode  voltage  pulses  respectively.  The  hollow  diode 

voltage  is  now  1.5MV,  just  as  in  the  experiment,  but  that  for  the  rod 
cathode  reaches  a  peak  value  of  about  1.7MV  which  is  0.7MV  too  great. 

This  was  adjusted  downward,  and  tested  as  reported  in  subsection  5.4  which  . 
follows.  The  final  diagnostic  shown  in  the  bottom  two  plots  in  the  figure 
is  the  B0  measurement  to  detect  the  presence  of  the  central  beam 
current.  Both  measurements  were  taken  at  a  radius  of  2.345  cm;  the  left 
at  Z  *  4.185  cm  and  the  right  at  Z  =  20.395  cm  (see  Figure  2).  They  both 
show  a  hollow  beam  measurement  of  about  -0.4  w/m2(  =  -4kG)  or  a  current  of 
about  47kA,  much  closer  to  the  experimentally  measured  value  than  the 
previous  case.  A  deflection  of  at  best  +0.025  w/m2(*  +0.25kG)  can  be 
detected  at  about  4  ns  indicating  a  central  beam  current  of  about  2.9kA. 
Since  the  peak  voltage  on  the  rod  cathode  is  so  exaggerated,  however,  this 
current  measurement  should  not  be  directly  related  to  the  actual 
experiment. 
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Figure  26.  Case  I  Time  Histories  for  (a)  Hollow  Cathode  Voltage,  (b)  Rod 
Cathode  Voltage,  (c)  B0  at  (34,24),  and  (d)  B0  at  (129,24) 

For  Dual  Beam  Case  I. 


5.4  Central  Beam  Only.  At  this  point  in  the  research  effort,  it  was 
deemed  necessary  to  simulate  the  CPA  model  device  in  the  absence  of  the 
hollow  beam.  Running  this  case  with  only  the  central  beam  present 

accomplished  two  major  objectives:  first,  the  applied  voltage  on  the  rod 
cathode  was  adjusted  downward  to  realize  the  correct  1.0  MV  of  measured 
voltage  and,  second,  a  benchmark  figure  for  the  full  current  of  the 

central  beam  was  obtained.  For  the  purposes  of  this  simulation,  the 

voltage  pulse  applied  to  the  rod  diode  was  changed  to  the  same  form  as 
that  which  had  been  applied  to  the  hollow  cathode  in  the  previous  runs. 
(No  voltage  is  applied  to  the  hollow  cathode  in  this  run.)  The  pulse 
turns  “on“  at  t  *  0  with  a  linear,  rapid  rise  to  a  plateau  value  of  0.615 
MV  in  0.0925  ns. 

The  formation  and  propagation  of  the  central  beam  for  these  operating 
conditions  is  documented  in  the  sample  electron  particle  plots  of  Figures 
27  and  28.  Some  degree  of  what  appears  to  be  striation  structure  is 

visible  in  the  downstream,  left  one- third  of  the  beam,  but  no  intensifi¬ 
cation  of  the  structure  or  perturbations  to  the  beam  propagation  is 
evident.  The  beam  seems  to  propagate  freely.  The  good  behavior  of  the 
beam  is  reinforced  by  the  classic  shape  of  the  axial  profile  shown  in 
Figure  29.  Measured  at  a  constant  radius  of  0.24  cm  (just  inside  the 
outer  radius  of  the  beam),  it  oscillates  only  mildly  about  zero  over 
almost  the  entire  length  of  the  beam  except  for  the  two  end  regions  where 
sharp  axial  acceleration  takes  place. 

The  validity  of  our  choice  for  the  applied  voltage  is  affirmed  by  the 
upper  left  time  history  plot  of  Figure  30.  A  near  perfect  1.0  MV  is 
measured  radially  across  the  rod  diode  gap.  Out  of  curiosity,  a  similar 
voltage  measurement  was  made  axially  down  the  length  of  the  tube  at  a 
radius  of  0.24  cm  between  the  rod  cathode  and  the  (floating)  hollow 
cathode.  This  is  shown  in  the  upper  right  plot  of  the  same  figure.  The 
transit  time  for  the  electromagnetic  wave  is  clearly  visible  in  the 
voltage  jump  at  about  one  nanosecond.  The  plateau  voltage  measured  in 
this  manner  is  only  slightly  less  than  1.0  MV.  To  gauge  the  effect  of  the 


Figure  27.  Central  Beam  at  t  =  500At  and  t  =  750At. 
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Figure  30.  Time  History  Plots  for  Radially  and  Axially  Measured  Voltages 
of  Rod  Cathode  (top  Left  and  Right)  and  Ez  Measured  at  tip 
of  Rod  Cathode  and  Hollow  Cathode  (bottom  Left  and  Right)  for 
Central  Beam. 
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central  beam  on  the  normal  electric  fields  at  the  tip  of  each  cathode 
which,  in  turn,  governs  the  emission  rate  there,  time  histories  of  E i 
measured  ‘at  the  tip  of  the  rod  cathode  (left)  and  hollow  cathode  (right) 
are  also  shown  in  the  bottom  two  plots  of  Figure  30.  A  quite  stable 
0.86  x  108V/m  is  measured  at  the  rod  cathode  while  only  an  insignificant 
1.0  x  106V/m  is  measured  at  the  hollow  cathode.  Clearly,  the  central  beam 
exerts  very  little  influence  on  the  emission  characteristics  of  the  hollow 
beam.  The  final  diagnostic  for  this  case  is  the  all-important  current 
measurement  benchmark  for  the  beam.  This  is  obtained  directly  from  the 
time  history  plots  of  Figure  31  which  track  B0  at  a  constant  radius  of 
1.41  cm  for  axial  distances  of  0.22  (top  left),  2.235  (top  right),  4.565 
(bottom  left),  and  19.955  cm  (bottom  right)  downstream  of  the  tip  of  the 
rod  cathode.  All  four  measurements  yield  a  plateau  B0  of  0.034w/m2 
(=  0.34  kG)  which  indicates  a  steady  central  beam  current  of  about  2.4 
kA.  This  is  in  excellent  agreement  with  experimental  findings  of  2-3  kA. 


5.5  Dual  Beam  Case  II.  There  now  remained  only  to  run  the  full,  dual 
beam  case  with  both  applied  voltages  properly  adjusted.  In  addition,  as 
in  the  previous  case,  the  rod  cathode  was  brought  linearly  to  full  voltage 
in  under  a  tenth  of  a  nanosecond  (after  the  hollow  beam  was  established) 
so  that  full  voltage  operation  of  the  central  beam  could  be  observed 
before  electrons  reflected  from  the  opposing  hollow  cathode  wall  could 
interfere  with  the  bulk  of  the  beam.  (Remember,  in  the  experiment,  the 
tube  was  5  meters  long  so  a  slow  pulse  rise  time  could  be  tolerated 
without  worrying  about  electrons  reflected  from  the  far  end  of  the  tube.) 
To  greatly  enhance  insight  into  the  detailed  physics  of  the  beams,  phase 
space  plots  were  also  added  to  the  simulation  diagnostics  package. 


As  in  the  previous  dual  beam  case,  the  hollow  beam  was  allowed  time  to 
equilibrate  before  a  voltage  was  applied  to  the  rod  diode.  Specifically, 
1.5  ns  elapsed  before  rod  cathode  "turn-on".  The  subsequent  development 
of  the  central  beam  is  glimpsed  at  t  =  2.45,  2.80,  3.15,  3.50,  3.85,  and 
4.20  ns  in  Figures  32  through  37  respectively.  The  left  half  of  each 
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Figure  31.  Time  History  Plots  for  8Q  Measured  at  R  =  1.41  cm  for 
Distances  of  0.22  (top  left),  2.235  (top  right),  4.565 
(bottom  left),  and  19.955  cm  (bottom  right)  Downstream  of 
Rod  Cathode  for  Central  Beam  Case. 
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Figure  32.  Real  Space  Plot  for  Both  Beams  and  Phase  Space  Plot  for 
Central  Central  Beam  for  Dual  Beam  Case  II  at  t  =  1750at. 


Central 


Figure  36.  Real  Space  Plot  for  Both  Beams  and  Phase  Space  Plot  for 
Central  Beam  for  Dual  Beam  Case  II  at  t  =  2750 At. 


Figure  37.  Real  Space  Plot  for  Both  Beams  and  Phase  Space  Plot  for 
Central  Beam  for  Dual  Beam  Case  II  at  t  =  3000At. 


figure  presents  a  sample  electron  position  plot  while  the  right  half  shows 
the  ?i  versus  Z  phase-space  plot  for  the  central  beam  electrons.  The 
key  result  to  be  gleaned  from  these  figures  Is  clear;  a  region  containing 
a  strongly  reflected  component  of  electrons  is  formed  almost  immediately 
In  the  central  beam  between  Z  *  23  and  Z  *  25  cm  (i.e.  right  off  the  tip 
of  the  rod  cathode).  This  indicates  the  presence  of  a  virtual  cathode 
there.  Such  an  entity  arises  due  to  a  local  build-up  of  electron 
space-charge  which,  in  turn,  probably  arose  from  the  previously  observed 
"reversed"  E- field  created  all  along  that  end  of  the  tube  due  to  the 
presence  of  the  outer  beam.  This  reflected  electron  structure  persists 
through  all  of  the  subsequent  time  figures.  The  corresponding  particle 
position  plots  also  show  evidence  of  enhanced  electron  number  density  in 
the  region  less  than  two  centimeters  downstream  from  the  rod  cathode  tip. 

In  all  but  the  first  of  this  set  of  figures,  are  visible  electrons 
reflected  from  the  opposing  hollow  cathode  wall  (which  is  at  a  lower 
potential).  These  electrons  clearly  wreak  havoc  with  the  extreme 
downstream  portion  of  the  central  beam  but  have  not  interfered  with  the 
important  upstream  physics  during  the  course  of  this  simulation.  The 
correspondence  with  the  experiment  is  therefore  complete. 

To  round  out  the  phase  space  picture  of  the  simulation,  two  further 
figures  are  provided.  Figure  38  plots  Pz  versus  Z  for  all  of  the 
electron  particles  in  the  hollow  beam  at  t  =  4.2  ns  into  the  simulation. 

It  shows  what  the  axial  phase  space  should  look  like  for  a  typical 

* 

foilless  diode  beam.  Contrasting  it  to  the  previous  plots  for  the  central 
beam  phase  space  highlights  how  severely  that  beam  has  been  disrupted  by 
the  presence  of  the  hollow  beam.  Figure  39,  on  the  other  hand,  shows  the 
central  beam  to  be  radially  "quieter"  than  the  hollow  beam.  It  plots  Pr 
versus  P i  side-by-side  for  both  beams.  Average  hollow  beam  electrons 
appear  to  have  a  significantly  greater  portion  of  their  kinetic  energy 
tied  to  radial  motion.  This  gives  rise  to  yet  another  interesting  point. 
Electrons  in  the  hollow  beam  with  finite  radial  velocities  will  interact 
with  the  25  kG  imposed  to  yield  a  finite  azimuthal  current  density. 
This  J9  (particularly  if  shank  emission  is  allowed)  could  result  in  a 
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Figure  38.  Phase  Space  Plot  (Pz  vs.  Z)  Beam  at 
t  *  3000&t  in  Dual  Beam  Case  II. 


Figure  39.  Plots  of  Pr  versus  Pz  for  Central  Beam  (left)  and 

Hollow  Beam  (right)  at  t  =  3000at  in  Dual  Beam  Case  II. 


diamagnetic  effect  Inside  the  hollow  beam.  This  could  thus  lower  the 
strength  of  the  axial  field  which  acts  upon  the  central  beam.  Distributed 
evenly  o\Jer  the  full  length  of  the  beams,  this  effect  is  probably  small 
for  the  present  case,  but  if  there  are  any  localized  axial  regions  of 
particularly  intense  J0,  then  the  collimation  of  the  central  beam  due  to 
the  applied  87  could  be  severely  disrupted.  The  present  simulations 
have  purposely  suppressed  any  azimuthal  motion  of  the  electrons  in  order 
to  limit  the  electromagnetics  to  the  first  mode  set,  (E^,  E^,  B0), 
thereby  streamlining  the  computational  costs.  Previous  diode  simulation 
experience  (see  References  17  and  18)  have  shown  it  to  be  very  safe  to 
Ignore  J0  in  parameter  ranges  such  as  these.  However,  if  the  decision 
Is  made  to  renew  the  CPA  experimental  effort,  new  simulations  should  be 
conducted  to  test  any  revised  and  improved  device  configuration  including 
J0  effects.  For  the  present  case,  the  observed  virtual  cathode  seems 
bad  enough  to  thwart  proper  operation.  Adding  JQ  effects  could  only 
make  things  worse  for  the  central  beam.  If  methods  are  found  to  mitigate 
this  virtual  cathode  effect,  then  a  check  on  J0  is  warranted. 

To  further  appreciate  the  weakening  effects  of  the  hollow  beam  on  the 
central  beam,  compare  the  axial  profile  of  Ej  just  inside  the  central 
beam's  outer  radius  for  this  run  shown  in  Figure  40  to  that  of  the 
previous  run  in  Figure  29  for  which  only  the  central  beam  was  present. 
The  Initial  accelerating,  field  at  the  far  right  has  been  cut  in  half. 
Furthermore,  It  is  followed  immediately  by  a  decelerating  region  (the  , 
virtual  cathode).  Then  the  field  oscillates  about  zero  until  hitting 
another  strong  decelerating  field  at  about  Z  =  9.0  cm  (probably  due  to  the 
electrons  reflecting  from  the  opposing  cathode  wall).  The  corresponding 
profile  for  the  hollow  beam  in  this  case  still  exhibits  the  typical 
foil  less  diode  beam  form  as  shown  in  Figure  41. 

The  voltage  characteristics  for  the  device  In  this  simulation  are 
shown  In  Figure  42.  The  upper  two  plots  display  the  time  histories  of  the 
voltage  pulses  on  the  hollow  cathode  (left)  and  rod  cathode  (right). 
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=  0.24  cm  (just  inside  central  beam)  at 
Case  II. 
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Figure  42.  Time  History  Plots  to  Voltage  on  Hollow  Cathode  (top 
left).  Voltage  on  Rod  Cathode  (top  right),  E7  near  Tip 
of  Hollow  Cathode  (bottom  left),  and  Ez  near  Tip  of  Rod 
Cathode  (bottom  right)  for  Dual  Beam  Case  II. 


The  plateau  hollow  cathode  voltage  is  extremely  close  to  the 
experimentally  measured  1.5  MV  while  that  for  the  rod  cathode  is  closer  to 
1.1  MV  cfempared  to  the  observed  1.0  MV.  The  lower  left  plot  traces  E^ 
as  a  function  of  time  at  the  hollow  cathode  tip  and  finds  it  to  be  quite 
steady  at  about  -0.9  x  108V/n,.  The  corresponding  bottom  right  plot  for 
E i  at  the  rod  cathode  tip  shows  a  brief  rise  to  0.5  x  108V/m  followed  by 
a  gentle  period  of  decay  to  0.42  x  108V/m. 

Figure  43  presents  the  key  diagnostic  to  this  entire  set  of  simulation 
runs,  namely  the  measurements  of  the  central  beam  axial  currents  in  the 

presence  of  the  hollow  beam.  The  locations  for  these  B0  measurements 

correspond  exactly  to  those  shown  in  Figure  31  for  the  central  beam 
wi thout  the  hollow  beam.  The  change  is  rather  dramatic.  The  sensors  at 
0.22,  2.235,  and  4.565  cm  from  the  rod  cathode  tip  all  show  brief  (»1  ns) 

rises  to  at  most  0.02  w/m2  followed  by  quick  declines  to  values  of  0.01 

w/m2  or  less.  The  bottom  right  plot  for  the  other  end  of  the  tube  detects 
no  significant,  consistent  axial  current  and,  in  fact,  seems  dominated  by 
low-level  displacement  current  noise,  probably  originating  from  the  hollow 
beam.  (The  presence  of  these  strong  dEjr/dt  oscillations  inside  the 

hollow  beam  argue  for  the  future  placement  of  the  B0  sensor  points  much 
closer  to  the  central  beam  to  help  sort  out  the  desired  signal  from  this 
background.  This  must  be  kept  in  mind  experimentally  too,  although  this 
computationally  observed  background  B0  is  probably  of  too  high  a 
frequency  even  to  be  detected  on  standard  experimental  hardware.)  The  , 
peak,  upstream  central  beam  current  that  these  measurements  indicate  is 
about  1.41  kA  which  quickly  decays  to  0.70  kA  or  less  in  approximately  a 
nanosecond.  This  current  is  down  by  over  a  factor  of  three  from  the  case 
without  the  hollow  beam.  It  explains  completely  why  no  central  beam 
propagation  was  detected  in  the  full  CPA  experiment  at  NRL. 


* 


77 


n/t\9*T 


m'jli  YEASI0H  JUNE  |<«3  OATE  10/76/09 
SIKJLAI 1(94  KftL  COLLECTIVE  PARTICLE  ACCELERATE 


MAGIC  VERSION  JUNE  !9fl3  O^TE  10/2A/P3 
SlHULflT itJN  NRL  COLLECTIVE  PARIJClE  ACCELEftATBR 


TIHE  HISTWT  PLOT 
e3  CCMPifNEMT 

*IQ-I  AT  CBCROIKATE  « 140. 151 


TI*F  HISTlWY  PLOT 
03  COMPONENT 

^jq*i  AT  CeOROINrtlE  1129.14) 


0.000  1.000  2.000  3.000  4.000  5.000 


TIME  (SEC)  ’  «1Q-* 


TIME  (SEC)  mIO*1 


KM1C  VERSI0N  JUNE  1983  OWE  10/26/63 
SIIU.RU8N  LSL  COLLECTIVE  WiKTIClE  RCCELtRRTOR 


vCKSliiN  jjtt  im  onrt  io//r,/m 

SHU.RII0I  NrtL  COLLECTIVE  PflRTICLE  RCCElERRTOR 


rri€  Hisrcrtr  n.ir 

83  CCNKOtCNT 

*10*'  »T  COORDINATE  ms.  Ml 


TIKE  I  SEC)  .10'’ 


*10-' 


rin-  mow.  nor 
03  UMfiWENI 
RT  CWRDimiE  139.  Ml 


TIME  (SEC)  *10'* 


Figure  43.  Time  History  Plots  for  B0  Measured  at  R  =  1.41  cm 

for  Distances  of  0.22  (top  left),  2.235  (top  right), 
4.565  (bottom  left),  and  19.955  cm  (bottom  right) 
Downstream  of  Rod  Cathode  for  Dual  Beam  Case  II. 
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6.  CONCLUSIONS 


As  a  Vesult  of  these  numerical  simulations,  it  can  be  concluded  that 
the  present  configuration  of  the  NRL  Collective  Particle  Accelerator 
causes  drastic  reduction  of  the  current  carrying  capacity  of  its  central 
beam  due  to  the  formation  of  an  effective,  virtual  cathode  immediately 
adjacent  to  the  physical  rod  cathode  tip  when  the  outer,  hollow  electron 
beam  is  present.  At  the  same  time,  however,  this  weakened  central  beam  is 
observed  to  propagate  stably  down  the  drift  tube.  It  seems  possible  that 
modest  redesign  may  be  able  to  raise  the  central  beam  current  to  a  level 
where  useful  CPA  experiments  could  be  conducted. 

Key  parameter  adjustments  which  might  lead  to  such  a  result  include  a 
lowering  of  hollow  cathode  voltage  (and  current)  and/or  a  raising  of  rod 
cathode  voltage.  Possible  changes  of  this  nature  must  first  be  evaluated 
by  the  experimentalists  in  light  of  their  effects  on  the  desired 
acceleration  mechanisms.  Proposed  adjustments  should  then  be  tested  again 
via  follow-on  simulations  to  confirm  their  positive  effects. 
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APPENDIX  A 


*  DISSPATCH:  A  Software  Package  to  Interface 

DISSPLA  Calls  with  NCAR  Graphics 

Mission  Research  Corporation's  two-dimensional,  electromagnetic, 
particle-in-cell  code,  MAGIC,  is  designed  to  operate  with  the  DISSPLA 
graphics  software  package  written  by  the  ISSCO  Graphics  Corporation.  The 
popularity  of  DISSPLA  nationwide  has  made  it  a  virtual  standard  in  the 
field  and  thus  the  most  logical  choice  for  use  in  connection  with  any 
transportable  code.  The  VAX/780  computer  available  to  Code  4780  at  the 
Naval  Research  Laboratory  does  not,  unfortunately,  have  the  DISSPLA 
library  available  to  it.  Instead,  it  uses  the  NCAR  Graphics  software 
package.  When  MAGIC  was  moved  to  the  NRL  VAX,  the  decision  was  made  to 
construct  a  software  library  which  translates  the  DISSPLA  subroutine  calls 
in  MAGIC  to  a  form  recognizable  to  the  NCAR  Graphics  library.  This 
approach  necessitated  the  construction  of  a  set  of  subroutines  whose  names 
corresponded  to  the  desired  DISSPLA  CALLS  and  which,  in  turn,  CALLed  the 
combination  of  NCAR  graphics  routines  which  most  closely  accomplished  the 
desired  DISSPLA  results.  The  set  of  new  subroutines,  written  by  MRC's 
James  McDonald,  forms  an  interface  package  which  has  been  dubbed, 
DISSPATCH.  Its  use  on  a  system  which  hosts  NCAR  graphics  eliminates,  as 
far  as  possible,  the  need  for  direct  modification  of  any  existing  code 
which  is  already  written  in  terms  of  DISSPLA.  The  price  the  user  must  pay  * 
for  this  convenience  is  the  failure  to  fully  utilize  the  real  power  of 
either  DISSPLA  or  NCAR.  Still,  the  plots  which  resulted  from  our  use  of 
DISSPATCH  were  very  pleasing  to  the  eye  and  looked  infinitely  better  than 
the  line-printer  graphics  which  we  would  have  otherwise  been  left  with. 

The  user  is  also  cautioned  that  the  present  version  of  DISSPATCH  only 
deals  with  a  subset  of  the  total  library  of  DISSPLA  calls.  Still,  that 
subset  is  large  enough  to  handle  the  likely  needs  of  most  users. 


Specifically,  the  MRC  MAGIC  code  required  NCAR  patches  for  the  following 

list  of  DISSPLA  calls: 

ANGLE 

AREA2D 

AXSPLT 

COMPRS 

CURVE 

DONEPL 

DOT 

ENOPL 

FRAME 

GRACE 

GRAF 

GRID 

HEIGHT 

MARKER 

MESSAG 

MRSCOD 

NOCHEK 

PAGE 

PHYSOR 

RESET 

SCLPIC 

THKFRM 

XMESS 

XNAME 

YNAME. 


These  25  separate  calls  were  subdivided  into  four  groups  according  to  the 
manner  in  which  they  were  handled  by  DISSPATCH.  The  first  group  was 
considered  "expendable."  These  were  deemed  to  be  not  important  enough  to 
warrant  the  trouble  of  writing  NCAR  counterparts.  This  group  included 
GRACE,  MARKER,  NOCHEK,  SCLPIC,  and  THKFRM.  This  was  coded  into  DISSPATCH 
using  the  set  of  dummy  subroutines  shown  in  Figure  A-l. 

The  second  group  consisted  of  all  DISSPLA  calls  which  have  NCAR 
counterparts  of  the  same  name  and  very  similar  function.  At  present, 
there  Is  only  one  member  in  this  group,  namely,  CURVE.  In  DISSPLA,  the 
call  to  CURVE  is  described  as* 

CALL  CURVE  ( XARAY ,  YARAY,  NPNTS,  IMARK) 


^  DISSPLA  Users  Manual,  Integrated  Software  Systems  Corp. ~ 

San  Diego,  California  (1981) 


c 

c 

c 

SUBROUTINE  GRACE (X) 
RETURN 
END 
C 

SUBROUTINE  MARKER! I) 
RETURN 
END 
C 

SUBROUTINE  NOCHEK 
RETURN 
END 
C 

SUBROUTINE  SClPIC(X) 
RETURN 
END 
C 

SUBROUTINE  THKFRM(X) 

RETURN 

END 


Figure  A-l.  The  "Expendable"  DISSPLA  Calls. 
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*  ^  V  jfu  n  mJl ilAx3  *1  itAiO.  O  nJ*  «J1  «wT.*0  !V , 


«* 


where 


*XARAY  =  name  of  array  containing  X  values 
YARAY  =  name  of  array  containing  Y  values 
NPNTS  =  number  of  points  to  be  plotted 
IMARK  =  frequency  of  marker  symbols. 


The  following  NCAR  call1,  differs  only  in  the  number  of  arguments  used: 


CALL  CURVE  (X,  Y.  N) 


where 

X  =  an  array  of  X  coordinates 
Y  *  an  array  of  Y  coordinates 
N  *  the  number  of  points  on  the  curve. 

Therefore,  the  DISSPLA  call  to  CURVE  can  be  used  as-is  with  the  IMARK 
argument  being  Ignored  by  NCAR  which  will  automatically  plot  every  point. 

The  third  "group"  also  contains  only  one  member  at  present.  This 
group  is  characterized  by  DISSPLA  calls  which  have  NCAR  counterparts  of 
the  same  name  but  different  functions.  The  subroutine,  FRAME,  fits  into 
this  group.  In  DISSPLA,  a  call  to  FRAME  draws  a  heavy  boarder  around  the 
subplot  area  while  In  NCAR  graphics  the  same  call  to  FRAME  advances  the 
plotter  one  frame  (similar  to  DISSPLA 's  ENDPL  call).  This  disparity  in 
functions  defied  any  simple  handling  with  DISSPATCH.  To  avoid  the 
redundancy  in  names,  a  subroutine  was  added  to  DISSPATCH  which  emulated 
DISSPLA' s  FRAME  with  NCAR  graphics  but  which  was  now  given  the  name, 
FRAMER.  Implementation  of  this  feature  with  MAGIC  required  the  changing 
of  every  occurence  of  "FRAME"  to  "FRAMER"  in  that  code.  A  listing  of  the 


t  The  System  Plot  Package,  NCAR  Technical  Note  NCAR-TN/162+1A, 

National  Center  for  Atmospheric  Research,  Boulder,  Colorado  (1981) 
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new  subroutine  is  provided  in  Figure  A-2.  It  references  a  graphics  COMMON 
block,  DISCOM  FOR,  which  is  shared  by  other  members  of  OISSPATCH  and  which 
is  shown*  in  Figure  A-3. 


The  final  group  of  subroutines  forms  the  bulk  of  the  DISSPATCH  package 
and  accurately  embodies  the  philosophy  of  emulating  DISSPLA  calls  with 
NCAR  graphics  "translations."  This  group  consists  of  the  following: 


ANGLE 

AREA20 

AXSPLT 

COMPRS 

OONEPL 

DOT 

ENDPL 

GRAF 

GRID 

HEIGHT 

MESSAG 

MRSCOD 

PAGE 

PHYSOR 

RESET 

XMESS 

XNAME 

YNAME . 


The  listings  for  these  subroutines  are  provided  in  Figure  A-4.  It  can  be 
seen  that  an  outside  subroutine,  PHYINT,  is  called  by  FRAMER,  GRID, 
HEIGHT,  and  MESSAG.  That,  utility  algorithm  is  included  in  DISSPATCH  to  . 
convert  DISSPLA  page  units  {inches  relative  to  the  physical  origin)  to 
NCAR ' s  plotter  address  units.  Its  listing  can  be  found  in  Figure  A-5. 

The  ultimate  question,  of  course,  is  "How  effective  has  the  DISSPATCH 
package  bsen  at  improving  the  graphics  output  by  MAGIC  at  NRL?"  The 
answer  is  best  given  by  presenting  samples  of  the  "before"  and  "after" 
plots.  Figure  A-6  is  a  sample  time  history  curve  as  generated  by  the 
line-printer,  graphics.  There  are  clearly  significant  problems  in  trying 
to  discern  the  "wiggles"  present  in  that  curve.  Fo-  comparison.  Figure 
A-7  shows  an  even  more  complicated,  longer  time  curve  using  DISSPATCH. 
Now,  even  the  smallest  "wiggle"  is  clearly  visable.  Another  key  MAGIC 


c 

SUBROUTINE  FRAMER 
INCLUDE  'OISCOM.FOR' 

C 

CALL  PHYINKO.,  0.,  MXO,  MYO) 
CALL  PHYINT{ AXX,  AXY,  MX,  MY) 
CALL  FRSTPT(MXO,  MYO) 

CALL  VECTOR{ MX,  MYO) 

CALL  VECTOR(MX,  MY) 

CALL  VECTOR(MXO,  MY) 

CALL  VECTOR! MXO.  MYO) 

RETURN 

END 


Figure  A-2.  The  NCAR  Replacement  for  OISSPLA's  FRAME 


I 


LOGICAL  PHYDEF 

COMMON  /DISPAR/  PAGEX,  PAGEY,  PFRACX,  PFRACY,  AXX,  AXY,  CHRWID, 

1  CHRHGT ,  ISIZ,  CHRANG,  IOR,  LXNAME(20) .IXNAME, 

2  LYNAME{2Q),  I Y NAME,  XOR,  XST,  XMX,  YOR,  YST,  YMX,  NDIVX, 

3  NOIVY,  PHX,  PHY,  PHYDEF,  LINETYPE 


Figure  A-3.  The  COMMON  Block  for  DISSPATCH. 
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o  o  u  o  o  o  o 


c 

SUBROUTINE  ANGLE(DEG) 
INCLUftE  'DISCOM.FOR' 

C 

CHRANG  =  DEG 
IOR  =  INT(DEG) 

RETURN 

END 


SUBROUTINE  AREAD2D(X,Y) 

INCLUDE  'DISCOM.FOR' 

AXX  =  X 
AXY  =  Y 

IF  (.  NOT  .PHYOEF)  CALL  PHYSOR(0.6  +  (PAGEX-AXX)/2.0,(PAGEY-AXY)/2.0) 

RETURN 

END 


SUBROUTINE  AXSPLT(AMIN,  AMAX,  AXMAX,  ORIG,  STEP.  AXIS) 
C 

*  DIMENSION  AMAN(4) ,  NDIVA{3) 

DATA  MANMAX  /  4  /  , 

1  AMAN  /  1.0,  2.0,  4.0,  5.0  /  , 

2  NDVMAX  /  3/  , 

3  NDIVA  /  4,  5,  6  / 

C  IGRT(Z)  *  INT(Z-AINT{Z)  +  1.0)  +  INT(Z)  -  1 
IF  (AMAX. GE. AMIN)  GO  TO  50 
STP  *  AMAX 
AMAX  *  AMIN 
AMIN  =  STP 
50  CONTINUE 

C 

BEST  =  0.0 
00  HD  =  1,  NDVMAX 
NDIV  *  NDIVA(ND) 

STP  =  (AMAX -AMIN) /NDIV 
LSTP  =  IGRT(AL0G10(ABS(STP) ) ) 

MSTP  =  0 


Figure  A-4.  The  DISSPATCH  DISSPLA  Emulators 


-* . 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS- 1963 -A 


100  MSTP  *  MSTP  +  1 

IF  (MSTP  .GT  .MANMAX)  THEM 
MSTP  -  1 

1  LSTP  *  LSTP  +  1 
ENO  IF 

STP  «  AMAM(MSTP)  *  10.0  **  LSTP 
RMIN  *  IGRT( AMIN/ STP)  *  STP 
RMAX  *  RMIN  +  STP  *  NDIV 
IF  (RMAX  .LT  .AMAX)  GO  TO  100 
GOODNESS  *  (AMAX -AMIN) /(RMAX-RMIN) 

IF  (GOODNESS  .GT  .BEST)  THEN 
BEST  »  GOOONESS 
ORIG  *  RMIN 
STEP  »  STP 
AXIS  *  FLOAT(NDIV) 

ENO  IF 

END  DO 
RETURN 
END 
C 
C 

c 

SUBROUTINE  COMPRS 
INCLUDE  'OISCOM.FOR' 

C 

•  CALL  SETK12,  12) 

CALL  PAGE(8.5,  11.0) 

CALL  HEIGHT( .14) 

CALL  ANGLE(O.O) 

RETURN 

END 

C 

C 

C 

SUBROUTINE  DONEPL 
CALL  FLUSH 
RETURN 
END 
C 
C 
C 

SUBROUTINE  OOT 
INCLUDE  'DISCOM.FOR* 

C 

CALL  OPTN( 'DPATTERN' ,  4681) 

RETURN 

END 

Figure  A-4  (continued). 
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VTiV? 


ooo  o  o  o 


SUBROUTINE  ENDPL(l) 
CALL  FRAME 
RETURN 
END 


SUBROUTINE  GRAFtXORIGD,  XSTPD,  XMAXD,  YORIGO,  YSTPD,  YMAXD) 
INCLUDE  'DISCOM.FOR' 

C 

DIMENSION  IFMTX(3) ,  IFMTY ( 3) 

C 

XORIG  =  XORIGO 
XSTP  »  XSTPO 
XMAX  »  XMAXD 
YORIG  «  YORIGO 
YSTP  a  YSTPD 
YMAX  *  YMAXD 

NDIVX  *  INT( (XMAX-XORIG)/XSTP  ♦  .5) 

NDIVY  «  INT((YMAX-YORIG)/YSTP  +  .5) 

XFAC  »  PFRACX/PAGEX 
YFAC  «  PFRACY/PAGEY 

.  CALL  SET ( PHY*XFAC , ( PHX  +  AXX)*XFAC.PHY*FAC.(PHY  ♦  AXY)  *  YFAC, 
1  XORIG,  XMAX,  YORIG,  YMAX,  1) 

IF  (IXNAME.EQ.O  .AND.  IYNAME.EQ.O)  RETURN 
CALL  FRAMER 

IXEX  -  I NT ( AL0G10( MAX( ABS ( XORIG) ,  kBS(XMAX)))) 

ENCODE! 10, 100,  IFMTX)  -IXEX 
100  FORMAT! ' ( ' ,  13,  ,PF6-3)1 > 

IYEX  a  INT(ALG010(MAX(ABS(Y0RIG),  ABS(YMAX)))) 

ENCODE (10, 100,  IFMTY)  -IYEX 

CALL  LABMOOt IFMTX,  IFMTY,  6,  6,  ISIZ,  IS1Z,  0,  2*ISIZ,  0) 

CALL  HALF AX( NDIVX,  1,  NOIVY,  1,  XORIG,  YORIG,  1.1) 

XL  «  XMESSUXNAME,  IXNAME,  0.,  0.) 

CALL  MESSAG(LXNAME,  IXNAME,  (AXX-XL)/Z.O,  -5.0*CHRHGT) 

IF  (IXEX.NE .0)  THEN 

CALL  MESSAG('*10\  3,  AXX-1.5*CHRWID,-5.0*CHRHGT) 

NC  -  1  +  I NT(ALG010( FLOAT (ABS (IXEX) )  +  .001)) 

IF  (NC.LT.l)NC  -  1 
IF( IXEX.LT.O)  NC  »  NC  ♦  1 
ENCODE (NC.  200,  EXP)  IXEX 
200  FORMAT ( I <NC>) 

CALL  HEIGHT ( CHRHGT /2 .0 ) 


Figure  A-4  (continued). 
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CALL  MESSAG(EXP,NC,AXX  +  3.0*CHRWID,  -9.0*CHRHGT) 
CALL  HE I GHT( CHRHGT*2 .0 ) 

END  IF 

YL  «  XMESS{LYNAME,  IYNAME,  0..  0.) 

CALL  ANGLE (90.0) 

CALL  MESSAG(LYNAME,  IYNAME.  -8.0*CHRWID,(AXY-YL)/2.0) 

CALL  RESET ( 'ANGLE* ) 

IF  (IYEX.NE.O)  THEN 

CALL  MESSAGE! '*10' ,  3.  -4.0*CHRWID,AXY  +  2.0*CHRHGT) 
NC  *  1  +  INT(AL0G10(FL0AT(ABS( I YEX) )  +  .001)) 

IF  (NC.LT.l)  NC  *  1 
IF  ( IYEX.LT.O)  NC  a  NC  +  1 
ENCOOEINC.  200.  EXP)  I YEX 
CALL  HEIGHT(CHRHGT/2.0) 

CALL  MESSAG(EXP,  NC.  -2.0*CHRWI0.  AXY  +  5.0*CHRHGT) 
CALL  HEIGHT(CHRHGT  *2.0) 

END  IF 
RETURN 
END 


SUBROUTINE  GRIO(IX.IY) 

INCLUDE  'DISCOM.FOR' 

’  DO  I  *  1,  NDIVX-1 

CALL  PHYINT(I*AXX/NDIVX,  -CHRHGT,  MX.  MY) 
CALL  FRSTPT(MX,  MY) 

CALL  PHYINT(I*AXX/NDIYX,  AXY.  MX.  MY) 

CALL  VECTOR! MX,  MY) 

END  DO 

DO  I  »  1,  NDIVY-1 

CALL  PHYINT(-CHRWID,  I*AXY/NDIVY,  MX,  MY) 
CALL  FRSTPTtMX,  MY.) 

CALL  PHYINT! AXX,  I*AXY/NDIYY,  MX,  MY) 

CALL  VECTOR! MX.  MY)‘ 

END  DO 
RETURN 
END 


SUBROUTINE  HEIGHT(HGT) 
INCLUDE  'DISCOM.FOR* 

CHRHGT  «  HGT 
CHRWID  «  HGT*6. 0/7.0 


Figure  A-4  (continued) 
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CALL  PHYINT(0.0,  0.0.  I.  J) 

CALL  PHYINT(CHRWID,  0.0.  12.  J) 

ISIZ  f  12-1 

RETURN 

END 


SUBROUTINE  MESSAG(ICHARS,  N.  X.  Y) 

INCLUDE  'OISCOM.FOR' 

DIMENSION  ICHARS(20) 

CHARACTERS*80  CHARS 

CONVERT  TO  UPPER  CASE. 

NM  *  1  +  (N-l)/4 

NRITE(CHARS,‘ (20A4) ' )  (ICHARS(I).  I  *  1.  NW) 

DO  I  »  1.  N 

K  *  I  CHAR  ( CHARS  ( I  :  D) 

IF  (K.GE.97  AND  K.LE.122)  X  *  K  -32 
CHARS( 1:1)  »  CHAR(K) 

END  DO 

CALL  PHYINT(X,  Y.  MX,  MY) 

CALL  PWRITIMX,  MY.  1REF( CHARS) .  N.  1S1Z.  10R,  -1) 
’  RETURN 
END 


SUBROUTINE  MRSCOO(TLENG,  NMRKSP,  RATR) 
INCLUDE  'OISCOM.FOR1 

CALL  OPTN( 'DPATTERN' .  255) 

RETURN 

END 


SUBROUTINE  PAGE(X,  Y) 
INaUOE  'OISCOM.FOR' 


Figure  A-4  (continued) 
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PAGEX  -  X 

PAGEY,  *  Y 

PMAX  «  MAX(X,  Y) 

PFRACX  «  PAGEX/PMAX 

PFRACY  «  PAGEY/PMAX 

RETURN 

END 


SUBROUTINE  PHYSOR(X,  Y) 
INCLUDE  'DISCOM.FOR' 

DATA  PHYDEF  /  .FALSE.  / 

PHX  *  X 
PHY  =  Y 

PHYDEF  «  .TRUE. 

RETURN 

END 


SUBROUTINE  RESET(X) 

-  INCLUDE  ‘DISCOM.FOR' 

DATA  DOT  /  '544F44'X  / 

IF ( X.EQ. ‘HEIG* )  CALL  HEIGHT! .14} 

IF(X.EQ.‘ANGL‘)  CALL  ANGLE (O.O) 

IF( X.EQ. DOT  .OR.  X.EQ. ‘DASH’ 

1  .OR.  X.EQ. ‘MRSC‘ )  CALL  OPTN( ‘DPATTERN* ,-l) 

RETURN 
ENO 


FUNCTION  XMES(LNESS,IMESS,DUM1,  DUM2) 
INCLUDE  'OISCOM.FOR' 

DIMENSION  LMESS(20) 

XMESS  *  CHRW1DMMESS 

RETURN 

END 


Figure  A-4  (continued) 
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SUBROUTINE  XNAME(LXNAM,  IXNAM) 
INCLUDE  'DISCOM.FOR' 

DIMENSION  LXNAM(20) 

I F ( IXNAM. NE.O)  THEN 

DO  I  «  1.  (IXNAM-D/4  +  1 
LXNAME(I)  «  LXNAM(I) 
ENO  00 
END  IF 

IXNAME  »  IXNAM 

RETURN 

END 


C 


C 


SUBROUTINE  YNAME(LYNAM,  1YNAM) 
INCLUDE  'DISCOM.FOR' 

DIMENSION  LYNAM(20) 

IF(IYNAM.NE.O)  THEN 

DP  I  «  1.  (IYNAM-D/4  +  1 
LYNAME(I)  -  LYNAM(I) 
END  00 

END  IF 

IYNAME  *  IYNAM 

RETURN 

END 


Figure  A-4  (concluded). 
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SUBROUTINE  PHYINT(X.Y.IX.IY) 

INCLUDE  'DISCOM.FOR' 

CALL  GETSIUEX.  IEY) 

IX  *  INT(2**IEX*(X  +  PHX ) / PAGEX*  PFRACX ) 
IY  *  INT(2**IEY*(Y  +  PHY)/PAGEY*PFRACY) 
RETURN 
END 


Figure  A-5.  The  PHYINT  Utility  Routine. 
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MAGIC  Vvrsicm  JUNE  1983  Dote.  7/29/83 
*  Simulation.  NRL  CHARGED  PARTICLE  ACCELERATUR 

TIME  HISTORY  PLOT 
Ji  COMPONENT 
AT  COORDINATE  <7,20 > 
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Figure  A-6.  Sample  Line-Printer  Curve 
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MAGIC  VERSION  JUNE  1983  DATE  9/11/83 
SIMULATION  NRL  C0LLECTIVE  PARTICLE  ACCELERATOR 


TIME  HISTORY  PLOT 
B3  C0MP0NENT 

*10-1  AT  COORDINATE  (34.24) 


Figure  A-7. 


Sample  DISSPATCH  Curve. 


diagnostic  consists  of  particle  position  plots.  The  best  approximation  a 
line  printer  can  make  Is  shown  In  Figure  A-8.  Much  useful  Information  Is 
given  in' the  plot  with  conducting  surfaces  labeled  with  C's  and  the 
symetry  boundary  labeled  with  S's.  In  addition,  the  spaclngs  between  the 
C*$  and  S's  correspond  to  the  actual,  variable  grid  spacing  used  In  the 
simulation  mesh.  However,  the  most  Important  part  of  the  plot  namely  the 
particle  positions  give  virtually  no  resolution  of  stream  fine  structure. 
This  Is  In  dramatic  contrast  to  the  DISSPATCH  particle  position  plot  of 
Figure  A-9. 

The  complete  source  listing  for  DISSPATCH  may  be  found  on  the  NRL, 
Building  A50  VAX/780  on  node  [BARKER. MCD]D I SSPLA. FOR.  The  corresponding 
object  code  Is  at  [BARKER.MCOlOISSPLA.OBJ.  It  Is  hoped  that  other  users 
will  exercise  the  package  and  make  suggestions  for  Its  Improvement  and/or 
expansion. 
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MAGIC  Version :  JUNE  1933  Date*  7/29/33 
Sinola t ion =  MRL  CHARGED  PARTICLE  ACCELERATOR 


TRAJECTORY  PLOT  OF  ELECTRONS  (  ISPE  -  1  ) 
AT  TIME  =  7.00E-11  SEC  FOR  i  TIrtE  STEPS 
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Figure  A-8.  Sample  Line  Printer  Particle  Plot. 
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MAGIC  VERSION  JUNE  1933  QRTE  6/17/83 
SIMULATION  NRL  COLLECTIVE  PARTICLE  ACCELERATOR 

TRAJECTORY  PLOT  OF  ELECTRONS  (  ISPE  =  1  ) 

AT  TIME  7.00E-10  SEC  FOR  1  TIME  STEPS 
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Figure  A-9.  Sample  DISSPATCH  Particle  Plot. 
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